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Abstract 

Alkaline-earth and ytterbium cold atomic gases make it possible to simulate SU(A^)-symmetric 
fermionic systems in a very controlled fashion. Such a high symmetry is expected to give rise to 
a variety of novel phenomena ranging from molecular Luttinger liquids to (symmetry-protected) 
topological phases. We review some of the phases that can be stabilized in a one dimensional 
lattice. The physics of this multi-component Fermi gas turns out to be much richer and more 
exotic than in the standard SU(2) case. For N > 2, the phase diagram is quite rich already in 
the case of the single-band model, including a molecular Luttinger liquid (with dominant su¬ 
perfluid instability in the A^-particle channel) for incommensurate hllings, as well as various 
Mott-insulating phases occurring at commensurate fillings. Particular attention will be paid to 
the cases with additional orbital degree of freedom (which is accessible experimentally either 
by taking into account two atomic states or by putting atoms in the p-band levels). We intro¬ 
duce two microscopic models which are relevant for these cases and discuss their symmetries 
and strong coupling limits. More intriguing phase diagrams are then presented including, for 
instance, symmetry protected topological phases characterized by non-trivial edge states. 

Keywords: alkaline-earth cold fermion, SU(A)-symmetry, Mott insulators, molecular 
superfluids, symmetry-protected topological phases 
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1. Introduction 

Symmetry lies at the heart of physics, since it is a powerful concept to classify conventional 
phases of matter (with broken symmetries) as well as phase transitions within the Ginzburg- 
Landau paradigm [1]. It has also been often used to idealize the actual physical situation, or 
simply to make analytical progress feasible as in the large-expansion. In high-energy physics, 
continuous symmetries based on the SU(iV) unitary group play fundamental roles in the stan¬ 
dard model of particle physics. For instance, an approximate SU(iV) symmetry, where N is the 
number of species of quarks, or flavors, underlies the description of hadrons. 
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In realistic condensed-matter experiments, such a high continuous symmetry is rarely re¬ 
alized [in stark contrast to SU(2) associated with, e.g., the spin-rotational symmetry] since it 
is not guaranteed by any fundamental principles of physics and usually requires some sort of 
fine-tuning of the parameters. There are, however, notable exceptions where relatively simple 
physical situations enable the system to fine-tune itself and realize a higher symmetry group. 
For instance, the interplay between spin and orbital degrees of freedom can lead to the realiza¬ 
tion of an SU(4) Kondo effect using semiconductor quantum dots [2] or an SU(4) symmetry in 
strongly correlated electrons with orbital degeneracy [3]. Another important example is electrons 
in graphene which have four flavors associated with the low-energy spin and valley degrees of 
freedom; the fractional quantum Hall effect in graphene is governed by the Coulomb interactions 
which are invariant under rotations of these four flavors and thus an SU(4) symmetry could also 
be realized [4, 5]. Though, weak atomic-range valley-dependent interactions might explicitly 
break this symmetry, it has been shown recently that an SO(5) symmetry arises in the i/ = 0 
quantum Hall regime of grapheme [6]. 

On top of these examples, ultracold fermions loaded into optical lattices might be ideal sys¬ 
tems to simulate strongly correlated electrons with a high symmetry, thanks to their great tunabil- 
ity (for reviews on recent development in many-body physics of cold atoms in optical lattices, 
see, e.g.. Refs. [7-9]). In principle, ultracold atomic gases with alkali atoms with hyperflne spin 
F" > 1/2 can explore the physics with SO(5) and SU(3) symmetries [10-15]. However, alkaline- 
earth atoms and related ones, like ytterbium atoms, have a peculiar energy spectrum associated 
with the two-valence outer electrons which make them the best candidates for experimental re¬ 
alizations of SU(iV) many-body physics [16-18]. The ground state (“p” state) is a long-lived 
singlet state ^S'o and moreover the spectrum contains a metastable triplet excited state (“e” state) 
^Pq. The g and e states have therefore zero electronic angular momentum, so that the nuclear 
spin / is almost decoupled from the electronic spin. From perturbation theory, the nuclear-spin- 
dependent variation of the scattering lengths is estimated to be smaller than ^ 10“® 10“^) 

for the g (the e) states [17]. As will be reviewed in Sec. 2.1.1, this decoupling of the electronic 
spin from the nuclear one in two-body collisions paves the way to the experimental realization of 
fermions with an SU(A^) symmetry where N = 21 + 1 (I being the nuclear spin) is the number 
of nuclear-spin states. 

The cooling of fermionic isotopes of these atoms below the quantum degeneracy has been 
achieved for strontium atoms ®^Sr with I = 9/2 [19, 20] (see Ref. [21] for a review) and for 
ytterbium atoms and respectively with I = 1/2 and 5/2 [22, 23] (see Ref. [24] 

for a review). These atoms enable the experimental exploration of the physics of fermions with 
an SU(A^) symmetry where N can be as large as 10. Several experiments have been done re¬ 
cently with these fermionic atoms loaded into optical lattices. An SU(6) Mott insulator has been 
explored with ^^®Yb atoms loaded into a three-dimensional (3D) optical lattice [25]. The experi¬ 
mental proof of the presence of the SU(A^) symmetry has been provided in Refs. [26] and [27] by 
using ®^Sr and ^^®Yb, respectively, in a 2D (®^Sr) and 3D (^^®Yb) optical lattices. Furthermore, 
the specific form of the interactions between the g and e states was determined in these works. 
Quite recently also, a coherent spin-orbital exchange interactions has been observed using ^^®Yb 
loaded into a 3D optical lattice [28]. 

All these remarkable experiments give us confidence that there will be many more upcom¬ 
ing breakthroughs in the near future. Motivated by recent experiments on ^^®Yb in the one¬ 
dimensional (ID) regime [29], in this review, we will be focusing only on the ID case in the 
presence of an optical lattice and will not discuss many other interesting results in two and three 
dimensions. Therefore, as a complement to our review, we refer the readers to Ref [18] for an re- 
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cent extensive review about the realization of SU(iV) symmetric fermionic gases and its physics 
in 2D or 3D. 

More specifically, we will hrst discuss the zero-temperature phases that can occur in the 
SU(A^) alkaline-earth fermions when only the atoms in the p-state are loaded in one-dimensional 
optical lattices. As will be seen in Sec. 2, the low-energy properties of alkaline-earth atoms 
in the 5 -state are described by the single-band SU(A^) Fermi-Hubbard model with the on-site 
interaction U that depends on the s-wave scattering length between two g atoms. When N > 
2 , rich physics emerges already in this simple setting, which is very different from that in the 
standard SU(2) case. At a hlling of one atom per site (1/A^-filling), which best avoids three-body 
losses, the repulsive S\J{N) Fermi-Hubbard model simplifies, for large repulsive interaction, 
to the S\J{N) antiferromagnetic Heisenberg spin chain with the fundamental representation at 
each site. The latter model, the so-called Sutherland model [30], has been quite extensively 
studied in the context of quantum magnetism [31-33] and the underlying gapless Mott-insulating 
phase can be directly investigated in ongoing experiments with ID ytterbium cold fermions. 
While low temperatures are difficult to achieve for fermionic gases, there is some advantage 
at working with several species (large N) since some short-range features might be observed at 
accessible temperatures due to some entropic effect, as was computed numerically in one [34, 35] 
and two dimensions [36]. The Mott insulating phase at 1/A^-filling has been already realized 
experimentally with ^^^Yb (N = 6 ) loaded in three-dimensional optical lattices [25]. 

In sharp contrast to the case of SU(2), the SU( A^) Fermi-Hubbard model displays/w/fy gapped 
Mott-insulating phases. However, these phases exhibits bond-ordering and spontaneously break 
translation symmetry. Therefore, the realization of some featureless gapped exotic phases of 
matter requires going beyond the single-band model. In this respect, it would be interesting 
to consider various two-band models which, on top of the S\J(N), have an additional ‘orbital’ 
degree of freedom. The origin of the two orbitals could be either the two atomic states g and e 
{g-e model), or the and py orbital states (p-band model) of the harmonic trap. 

In the two-orbital cases, the interplay between the orbital and SU(A^) nuclear-spin degrees of 
freedom will be shown to give rise to several interesting phases, including symmetry-protected 
topological (SPT) phases [37, 38]. The latter refer to non-degenerate fully gapped phases which 
do not break any symmetry and defy the characterization with local order parameters. Since 
any gapful phases in one dimension necessarily have short-range entanglement [39] and can be 
reduced to a trivial product state by a series of local unitary transformations, the presence of a 
symmetry (which restricts the type of possible local unitary) is necessary to protect the properties 
of that ID topological phase, in particular the existence of non-trivial edge states [38, 39]. 

The outline of this review is as follows. In Sec. 2, we will provide a minimal introduc¬ 
tion to SU(A^) symmetry in alkaline-earth systems and then describe the known results for the 
single-band Hubbard model in ID. We consider both the cases with commensurate hllings, where 
various Mott phases occur, and those with incommensurate hllings where low-energy properties 
can be described using the Luttinger-liquid theory [40, 41]. 

In Sec. 3, we will move on to the case of two-‘orbital’ fermions. After introducing two 
relevant microscopic models (g-e model and p-band model), which are directly relevant to ex¬ 
periments [26, 27], we will discuss their strong-coupling limits, i.e., the S\J(N) spin chains in 
various representations, for which several results are known or conjectured. In particular, we 
will sketch the known classihcation of the so-called SU(W) SPT phases allowing us to under¬ 
stand the formation of topological states in our SU(A^) fermion system as generalizations of the 
well-known Haldane phases in spin systems. Finally, we will present some typical numerical 
results to provide complete phase diagrams as well as clear signatures of the SPT phases through 
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the direct measurements of nontrivial edge states. Conclusions will be given in Sec. 4. 


2. Single-band Fermi-Hubbard model 

2.1. Definition of the model 

2.1.1. SU(N) symmetry in cold-atom systems 

Like the SU(3)-symmetry in the quantum chromodynamics for the strong interactions, in 
high-energy physics, the existence of the multiplet of N particles and the SU(iV)-symmetry that 
governs the dynamics of these particles have nothing to do with spins and are built in the the¬ 
ory from the outset. In condensed-matter physics, on the other hand, such internal symmetries 
independent of spin may originate from, e.g., sublattice symmetry and normally fine-tuning is 
necessary to realize high symmetries as discussed in the introduction. In order to understand 
how SU(iV) symmetry arises among the multiplet which originates from the angular momentum 
without any fine-tuning, consider two spin-F" fermionic atoms interacting with a contact interac¬ 
tion. Due to the rotational invariance, the two-body interaction depends only on the total spin /. 
Among all the possible values / = 0,1, 2,..., 2F, only the {2F -f 1)/2 antisymmetric combi¬ 
nations / = 0, 2,..., 2F — 1 (F is assumed to be half-odd integer) are allowed by fermionic 
statistics. Therefore, the most general form of rotationally [i.e., SU(2)] invariant interaction 
between two spin-F objects may be written as [42] ' 


2F-1 



( 1 ) 


/= 0 , 2 ,... 


where Pf denotes the projection operator onto the total spin-/ sector and the coupling constants 
gf are determined by the corresponding s-wave scattering lengths a/ as (M being the mass of 
the atoms in question) 


In the ground state of for instance, I = F = 5/2 and we have, in principle, three 

independent coupling constants go, g 2 , and g^. 

In the second-quantized formulation, the above projection operators can be simply given as 


/ 


/ 



( 3 ) 


m=-f 


m=-f 


with the ‘pairing operators’ defined by [42] 


F 



(4) 


F 



^ A similar analysis has been done in Ref. [43] for the F = 1 bosons. 
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Using the definition of the pairing operators [i.e., creation/annihilation operators for a pair with 
a particular value of the hyperfine spin (/, m)], the interaction may be rewritten as [42, 44] 

2F-1 

2 SfPf 

/=0,2,... 

1 2F-1 IFF (5) 

= 2 ^ ^ ^ ^ {F, ai;F, / 3 i|/, m){f, m\F, 02; F, P2) 

/— 0 , 2 ,... m— — f 0 . 1 ,01 —— F O2i02 — — F 

X 4i(l')4iW%WCa2(l’) ■ 

When the coupling constants (or, the scattering lengths) do not depend on / (i.e., gf = g), the 
above simplifies to the following expression 

2F-1 2F-1 / 

2^ H ^f= 2^ £ ^/.m(0^/.m(r) 

/=0,2.... /=0,2,...m=-/ 

1 ^ ^ 

= E E 

oi,0i — — F 02,02——F 

{ 2F-1 f ^ 

E E {F,0!i;F,Pi\f,m){f,m\F,a2;F,l32) > (6) 

f=0,2,... m=-f J 

'■-V-" 

— ,02 ,132 

X 4i(l')4iWc/52(r)Ca2(l’) 

1 ^ 1 
= 2 ^ E 4 (i’) 4 Wc/ 3 (i')ca(r) = -5 : n(r)n(r) : , 

0,0—— F 

where n(r) = cj,(r)ca(r) is the density operator which is invariant under SU(7V) symmetry 
with N = 2F + 1. Therefore, when gf do not depend on /, i.e., when the fine-tuning of the 
scattering lengths oq = 02 = ■ • • = a 2 F-i occurs, the original SU(2)-symmetry of the atom- 
atom interaction gets enlarged to SU(iV) with N = 2F -f 1^. 

A natural question to ask then is when and in which system this fine tuning happens? The col¬ 
lision among neutral atoms are mainly governed by short-range van der Waals interaction that de¬ 
pends only on the electronic wave functions of the two colliding atoms. The nuclear-spin degrees 
of freedom can participate in the collisional processes only through the hyperfine interaction. In 
alkali cold atoms with hyperfine F = 3/2 spin, an enlarged continuous Sp(4) ~ SO(5) symme¬ 
try arises without fine-tuning^, i.e., with two independent scattering lengths oq, a 2 - [10, 13] The 


^Even when the fine-tuning is incomplete, we may have other extended symmetries. See Ref. [45] for the discussion 
of the phase structure of these cases. 

^In order to understand this, we first identify the two quartets 4 of the two colliding F = 3/2 atoms with two 4- 
dimensional spinor representations (4so(5)) of S0(5). Then, the states of the two atoms may be decomposed, in terms of 
S0(5), as 

{F = 3/2) (gi (F = 3/2) 4so(5) ® 4so(5) — lso(5) © 5so(5) © 10so(5) , 

among which the last one (symmetric 10so(5)) must be discarded by the Fermi statistics (for the s-wave scattering). The 
remaining two (lso( 5 ) 5so(5)) may be identified with the singlet (/ = 0) and the quintet (/ = 2) of SU(2) [see 

Eq. (1)]. Therefore, the scattering of two fermionic F = 3/2 atoms can be rephrased in terms of S0(5). 
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physical properties of this model in one dimension have been extensively investigated over the 
years by means of various analytical and numerical approaches [12, 45-56]. The fine-tuning of 
the scattering lengths oq = a 2 to achieve a higher SU(4) symmetry is not easy to realize exper¬ 
imentally with these alkali atoms. In contrast, if the hyperfine interaction is quenched for some 
reasons, the scattering processes are independent of the nuclear spin /: 5o = 52 = • • • = 52 /- 1 - 
As the total electron angular momentum vanishes (and so does the hyperfine interaction) in the 
ground state of alkaline-earth and Yb atoms, these atoms are the best candidates of systems 
that realize SU(iV) symmetry without any fine-tuning [16-18]. 


2.1.2. SU(N) Hubbard model 

In order to derive a lattice Hamiltonian that describes low-energy physics of an A^-component 
Fermi gas with the SU(7V) symmetry moving on ID optical lattices, let us begin with the single¬ 
particle problem: 

no = + ^per(^)| = n//{z) , (7) 

where V.pe^{z) is a periodic potential that introduces a lattice structure along the chain (i.e. z) 
direction. Also we have assumed that the confining potential in the other directions (i.e., x and 
y) is so strong that we are able to neglect the motion in the xy-direction and set x = y = 0. If the 
chain is translationally invariant in the z-direction, the single-particle state is given by the Bloch 
function (z) (n being the band index) satisfying: 

n/i{z)ip^^^{z) = . ( 8 ) 

To derive an effective Hubbard-like Hamiltonian [57, 58], it is convenient to move from the 
Bloch function (z) to the Wannier function defined by 

vAcell ^ 


(the index R labels the center of the Wannier function and A];eii is the number of unit cells in the 
z-direction) and introduce the fermion operators in the Wannier basis 

Ca(^) = > 4(^) = ^ 

R n R n (10) 

(a = l,...,iV) . 


We use these operators to rewrite the SU(Af)-invariant two-body interaction (6) as 


Y 

a,/3=-F' 


dz cl (z)ci (z)Cf3 (z)Ca (z) 


_ T//ni.n2;Ji3,"4 \ (n2)f (ns) ( 714 ) 

2 2^ 2^ (Ri,R2-,R3,RA>^a.Ri^l3,R2''P,R3^a,R4 ' 

a,f3=l {Ri,ni} 


( 11 ) 


where 


dz' 


,,("1 

^Ri 


iz)w^l^>iz)wZ^>iz) 


(ns). 


,,("4)/ 


-\/-f‘ni,n2;n2,n4 \ ^ 


6 


( 12 ) 




and we have re-labelled a = —F,... ,F ^ a = 1,... ,N (N = 2F + 1)"^. As usual [57, 58], 
we keep only the most relevant band (denoted by nx = = ri'i = = Uq) and the on-site 

term Ri = R 2 = R 3 = R 4 = Ri to obtain the so-called Hubbard interaction: 


1 ^ 

^Y(no,no-,no,no-, „t t . 


N 


(13) 




a.,^—1 i 5 

In the above, we have introduced a short-hand notations 

i Ri 


(14) 


and rii = '^1 ^Ca, t is the density operator on site i. 

When second-quantized, the single-particle part reads as 

AT 

^ dzci{z)n//{z)ca(z) 

= { hzw^^^'>*{z)Hi,{z)w^^^\z)\ 

a=l{R,,ni} ^ 

This may be further rewritten using the hopping amplitudes: 


(ni)t (n2) 
^Ol,Ri ^Q:,i?2 ’ 


(15) 


J dzw^R^^*{z)'H//{z)w^^^\z) 

kz k'^ 

\ = - R2) , 



and we finally obtain the kinetic term 

-E E E 

Q!=l nGbands Ri ,i?2 


(16) 


(17) 


Retaining only the terms with n = uq and | i?i — i? 2 1 = 1 in the above and combining it with 
Eq. (13), we arrive at the S\J{N) generalization of the famous Fermi-Hubbard model: 

N ^ 

RsV(N) = —t EE( *-|-l “1“ H.C.^ + ^ ^^ Tlj ijli 1) . (18) 

i Q!=l i 


fact, N can take any number N < 2F -h 1 as one can load only a subset of the multiplet by optical pumping. 
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The model (18) is invariant under the global U(l) symmetry: Caj i—> e^^Caj, which implies 
the conservation of the total number of atoms. In the rest of this paper, we frequently use the 
terminology ‘charge’ to denote the degree of freedom associated with this symmetry, although 
we are dealing with neutral atoms without electric charge. All the parameters in the model (18) 
are independent from the nuclear-spin states (a = 1,..., N) and an extended SU(IV) symmetry 
arises: Ca, i i— i, U being an S\J(N) matrix. The actual continuous symmetry group 
of the Hamiltonian (18) is then U(iV) = U(l) x SU(A^) but the model "HsutAT) is often called 
the SU(IV) Fermi-Hubbard model to put the emphasis on its non-trivial SU(iV) hyperfine-spin 
rotational invariance. 

2.2. Sutherland model and its low-energy physics 

The model (18) describes alkaline-earth atoms in the g state (i.e., ground state ^S'o) loaded 
into the lowest band of the optical lattice. The interaction parameter U is directly related to 
the s-wave scattering length associated with the collision between two atoms in the g state [see 
Eqs. (12) and (13)]. When N = 2, the model (18) is the usual SU(2) Hubbard chain which 
is exactly solvable by means of the Bethe ansatz [59]. The physical properties of the model 
have been discussed in great detail over the years and are reviewed in the book [60]. However, 
for N > 2, the Hamiltonian (18) is not integrable for arbitrary U and filling n. Although it is 
possible to formally generalize the Lieb-Wu Bethe ansatz equation [59] to fermions with internal 
SU(A^) symmetry, it is believed that the corresponding model describes a non-local variant of 
the SU(A^) Hubbard model [61]. In the absence of a lattice, the model is again integrable and its 
properties have been described in a recent review [62]. 

The situation becomes much simpler in the limit of large repulsive U for a filling n = 1 /N 
with one atom per site which best avoids the three-body losses. In that case, the model (18) 
reduces to the S\J(N) Heisenberg antiferromagnetic spin chain with the SU(A^) fundamental 
representation (represented by the Young diagram □; for a pedagogical explanation of the repre¬ 
sentation theory of SU(A^) and the Young diagrams, see, e.g.. Ref. [63]) on each site {Sutherland 
model [30]): 

H = (19) 

i 

where J = 2t^ jU is the antiferromagnetic spin exchange and is the operator which 

permutes the Sl]{N) hyperfine states on the sites i and For the fundamental representation 
(□), Pi. i+i is compactly written, in terms of the SU(iV) generators that are normalized to be 
TrSfSP = 6^^, as 

N^-l 

= n+T. ^20) 

A=1 

[SU(iV) generalization of the Dirac identity]. It is known that the ‘spin’ model (19) well describes 
the low-energy sector of the original Hubbard model (18) for U/t > 12 [64]. 

In contrast to the original fermionic model (18), the large-C/ effective Hamiltonian (19) can 
be solved exactly by the Bethe ansatz [30]. The low-energy spectrum is gapless with N — 1 
relativistic modes with the same velocity = 2ttJ/N. The critical theory has been identified 
by Affleck [31, 32] as described by the level-1 SU(W) Wess-Zumino-Witten [SU(7V)i WZW] 
conformal field theory (CFT) [65-67]. This CFT has a central charge c = N — 1 and the low- 
temperature specific heat (per volume) scales as [68, 69]: C{T) ~ l)T/(6hJ) (with 
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/cb denoting the Boltzmann constant). The latter result has been conhrmed by the thermody¬ 
namic Bethe ansatz [70] (see Ref. [33] for the determination of c by finite-size corrections) and 
quantum Monte-Carlo calculations (QMC) [35, 71]. It is known [32] that the SU(7V)i WZW 
CFT corresponds to the stable hxed point of generic ID gapless systems with S\J{N) symmetry 
[as the SU(2)i WZW CFT describes generic gapless SU(2)-invariant spin chains]. 

The gapless behavior in the SU(7V) Mott-insulating phase with one atom per site manifests 
itself in the spin-spin correlation functions which exhibit a universal power-law decay in the 
long-distance limit at zero temperature [31, 32]: 



( 21 ) 


where r is the Euclidean time, gq is the lattice spacing, and A,B = 1,... — 1 are the 

components of the SU(A^) spin operator. The logarithmic corrections in Eq. (21) stem from the 
existence of a marginal operator in the low-energy effective field theory describing the model 
(19) [72, 73]. When this marginal operator becomes marginally relevant by adding, e.g., fur¬ 
ther neighbor interactions, the system may spontaneously develop the A^-merization [i.e., the 
formation of iV-site clusters in SU(A^)-singlet] and enter gapped phases with broken translation 
symmetry [74—76]. 

The thermodynamics properties of the Sutherland model (19) have been investigated numeri¬ 
cally in Refs. [34, 35]. In particular, it has been shown that characteristic short-range correlations 
develop at low temperature as a precursor of the algebraic correlations (21) in the ground state. 
The hrst sign of short-range order appears at an entropy per particle which increases with TV, 
leading to observable qualitative effects in ongoing ultracold atom experiments with alkaline- 
earth fermions. 

2.3. Phase structure 

Since the SU(iV) Eermi-Hubbard model (18) with N > 2 is, in general, not exactly solv¬ 
able, one has to resort to approximate but powerful techniques available in one dimension to 
map out its zero-temperature phase diagram: field-theoretical, strong-coupling, and numerical 
approaches. In this section, we sketch the phase diagram of the model (18) for both incommen¬ 
surate (Sec. 2.3.1) and commensurate fillings (Secs. 2.3. 2-2.3.4). The main results of Sec. 2.3 
are summarized in Table 1 . 

2.3.1. Incommensurate fillings 

We first consider the SU(iV) Eermi-Hubbard model (18) for incommensurate fillings. In this 
respect, there are no umklapp processes which open a gap for the charge degrees of freedom [or, 
the U(l) sector corresponding to the continuous symmetry of Eq. (18)]. A metallic state is then 
formed whose nature strongly depends on the sign of the coupling constant U. 

When U > 0, all modes are gapless and a metallic Luttinger-liquid phase emerges [40, 41]. 
The hallmark of this TV-component Luttinger liquid phase is the 2kfi oscillations in the density- 
density and the spin-spin correlation functions decaying with non-universal power-law exponents 
[32, 64, 77]. Eor instance, the leading asymptotics of the SU(TV) spin-spin correlation functions 
(21) reads now in the metallic phase as: 


(5^(r,a:)5^(0,0)) ~ 


cos {2kfx) 



( 22 ) 





where K and Vc are respectively the Luttinger parameter and the characteristic velocity of the 
charge excitation which depend on the interaction and density [40, 41]. In particular, K de¬ 
termines the singularity in the momentum distribution around the Fermi point k^: n{k) ~ 
n{k-p) + const, sign (fc — fcp) \k — A:f|“ with a = (1 — K)"^/{2NK) [77]. This power-law 
singularity at the Fermi level is inherent in Luttinger liquids unlike in the standard Fermi liq¬ 
uid [40, 41]. Similarly, the single-particle density of states also has an anomalous power-law 
behavior for any finite value of N: p{uj) ^ |u;|“. 

On the other hand, the physics is very different in the attractive case {U < 0); there is a spin 
gap for the SU(iV) degrees of freedom and the only remaining gapless mode is the charge one. 
The metallic phase is then characterized by a CFT with the central charge c = 1 reflecting a single 
bosonic gapless mode. The resulting spin-gap phase is the Luther-Emery phase [40, 41, 78], 
which turns out to be very exotic when > 2 at sufficiently low density [50, 79]. In the 
N = 2 case, the Luther-Emery phase describes the competition between a charge-density wave 
(CDW) instability and a superconducting one [40, 41]. In the attractive SU(2) Hubbard model, 
the leading instability is the superconducting one. When > 2, on the other hand, the SU(A^) 
symmetry plays an important role in one dimension by preventing any pairing between fermions: 
there is no way to form an SU(A^) singlet with only two fermions. When N > 2, the usual 
pairing instability is then completely suppressed with exponential-decaying correlation functions 
in stark contrast to the N — 2 case. The only possible gapless fluctuation corresponding to the 
superfluid instability is a molecular one where N fermions form S\J(N) singlet (an analog of 
baryons in high-energy physics): J\4i = c\ jcj ^... cjy j, i.e, trionic (N = 3) and quartetting 
(N = 4) superfluid instabilities. The Luther-Emery phase of the model (18) with (7 < 0 is then 
governed by the competition between the instability toward this molecular superfluid (MS) and 
the one toward CDW. Among these competing order parameters, the one that exhibits the slowest 
power-law-decaying correlations at zero temperature corresponds to the leading instability. The 
equal-time correlation functions of these order parameters have been determined by means of 
bosonization [12, 49, 79] as: 

{nim+x) ~ cos {2kpx) 

I f for N even (23) 

* |sin (/cpa:) for A^ odd. 

Either CDW or MS instability thus dominates depending on the value of the Luttinger parameter 
AT; dominant MS instability requires K > N/y/3 (K > N/2) when N is odd (even). At issue is 
the value of the the Luttinger parameter K. Its expression as a function of U and the density n 
has been numerically determined in Ref. [79] and is reproduced in Eig. 1 for A^ = 3. In the low- 
density regime, the MS phase, characterized by the bound states made of N fermions, exists for 
a wide range of attractive U (the shaded region in Eig. 1). The latter phase might also be viewed 
as a ‘molecular’ Luttinger liquid with the molecules of N atoms and, as has been mentioned 
above, is characterized by the suppression of the usual Cooper pairs. In the high-density regime 
(the lower part of Eig. 1), on the other hand, the dominant instability for [/ < 0 is the standard 
CDW. 

2.3.2. Mott transition 

Eor commensurate fillings, umklapp processes might become strongly relevant perturbations 
and open a spectral gap for the charge degrees of freedom. The nature of the resulting Mott- 
insulating phase depends crucially on the Ailing n [80]. In the case of one atom per site (i.e., 
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Figure 1: (Color online) The Luttinger parameter K as function of U and the density n for the SU(3) Fermi-Hubbard 
model with attractive interaction (U < 0). Grey region marks the onset of the MS phase. The n = 3l2 line corresponds 
to the half-filled case where a fully gapped Mott-insulating phase occurs as discussed in Sec. 2.3.3. From Ref. 179]. 


n = 1/N filling), the Mott-insulating phase is gapless in the large-{7 limit and is described by 
the Sutherland model (19) and its physical properties were discussed already in Sec. 2.2. In the 
N — 2 case, it is well-known from the exact Bethe-ansatz solution that there is no (hnite-[/) 
Mott transition and a charge gap opens exponentially as soon as positive U is switched on [59]. 
A gapless Mott-insulating phase with central charge c = 1, described by the SU(2) Heisenberg 
spin model, emerges thus for all U > 0. 

Again, the situation might be very different when N > 2 as advocated in Ref. [77]. In 
the weak-coupling (U <C t) limit, the umklapp operator appearing in the low-energy effective 
Hamiltonian has the scaling dimension A — NK [77], which is relevant when the Luttinger 
parameter is sufficiently small; K < 2/N. When N = 2, this is always the case since K < 1 
as soon as a repulsive interaction is switched on [40, 41]. In stark contrast, the umklapp operator 
is always irrelevant when N > 2 for sufficiently small U and a Mott-transition at hnite U is 
thus expected toward the large-(7 insulating phase with the gapless S\J{N) sector described by 
the SU(iV) Sutherland model (19) [77]. This problem has been investigated numerically by 
various different techniques [64, 77, 81, 82]. The determination of the position of the Mott- 
transition turned out to be a difficult numerical problem due to the smallness of the charge gap. 
The QMC simulations of Ref. [77] found a critical value Uc — 2.2t (N = 3) and Uc — 2.8t 
(N = 4). Density matrix renormalization group (DMRG) [83] calculations done in Ref. [64] 
reported smaller values: (7c ~ l.lf (N = 3) and Uc — 2.It (N = 4) while those in Ref. [82] 
concluded the same value Uc = 0 for all N > 2. The physical properties of the metallic phase 
when U < 64 are described by the multi-component Luttinger liquid with N gapless channels 
(see Sec. 2.3.1). Above the critical value U > Uc, there is a charge gap and we are left with 
N —1 gapless spin modes whose physical properties are governed by the Sutherland model (19). 
The Mott-transition is argued [77] to be of the Berezinskii-Kosterlitz-Thouless universality class 
[84, 85]. 
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2.3.3. Half-filled case 

We now consider the half-filled case with fcp = 'Kl2aQ (N/2 fermions per site; n = 1/2). 
For = 2, it corresponds to the situation where we have one atom per site and reduces to the 
case considered already in Sec. 2.2; the physics for [/ > 0 is governed by the Heisenberg model 
with a gapless c = 1 behavior [corresponding to the level-1 SU(2) WZW CFT]. On the attractive 
side U < 0, we can apply a transformation (Shiba transformation [86, 87]; See section 2.2.4 of 
Ref. [60] for a detailed discussion of the Shiba transformation), that interchanges spin and charge 
while flipping the sign of U, to show that now a gap opens in the spin sector while the charge 
sector remains gapless (Luther-Emery liquid). 

In contrast, when N > 2, all degrees of freedom are fully gapped for any values of U 
(whether positive or negative) due to the absence of spin-charge separation in the low-energy 
limit [45, 82]. The resulting fully gapped Mott-insulating phase is two-fold degenerate as the 
result of the spontaneous breakdown of the one-site translation symmetry [45]. The physical 
nature of the Mott-insulating phase depends crucially on the sign of U. 

In the attractive case {U < 0), long-range ordering of period-2 CDW emerges [45, 88, 89]. 
In the strong-coupling region (\U\ ^ t), the picture of the CDW formation is simple; SU(7V)- 
singlet molecules of N atoms (A^-mers) are formed first and then the preformed A^-mers organize 
themselves into period-2 crystalline structures in such a way that they optimize the 0(t'^/\U\) 
repulsive interaction generated by virtual hopping [88, 89]. To illustrate the crystalline pattern 
of A^-mers, we show, in the left panel of Fig. 2, the spatial profile of physical quantities (lo¬ 
cal fermion density and kinetic-energy density) obtained by DMRG for the half-filled SU(4) 
Hubbard model at U/t = —8. The local fermion density clearly shows period-two oscillation 
indicative of CDW of tetramers (note that the maxima of the density is close to 4), while the 
kinetic energy does not exhibit any special feature. A typical wave function of the CDW phase 
is shown in Fig. 3(a). 



Figure 2: (Color online) Local density and kinetic bond energy for the fermionic SU(4) model with two particles per 
site at (7/f = —8 (left, data are only shown in the bulk) and U/t = 8 (right) on a L = 64 chain computed by DMRG. 
Data strongly indicate CDW and SP phase respectively. 


When U > 0, on the other hand, it has been shown that a gapful dimerized [or spin-Peierls 
(SP)] phase with bond ordering [see Fig. 3(b) for an intuitive picture of the ground state] appears 
upon switching on the weak repulsive interaction [45, 82]. In this respect, the situation is very 
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different from the N = 2 case where, at half-filling (one atom per site), a Mott-insulating phase 
with gapless spin excitations is stabilized when ?7 > 0. It has been shown numerically (by means 
of QMC and DMRG) for TV = 4 that there is adiabatic continuity from weak to strong coupling 
and that the SP phase occurs for all 17 >0 [45, 90]. In the large-t/ limit, the existence of this SP 
phase can be simply understood from the fact that the half-filled SU(4) Fermi-Hubbard model 
reduces to the SU(4) Heisenberg spin chain in the antisymmetric self-conjugate representation 
(Q) of SU(4); the latter model is known to display, at zero-temperature, a dimerized phase with 
two-fold ground-state degeneracy [91-93]. In the right panel of Fig. 2, we also present the plot 
of the local fermion density and the kinetic-energy density for the SP phase obtained by direct 
DMRG simulations of the SU(4) Hubbard model at U/t = 8. As can be clearly seen, there is 
a period-two oscillation in the profile of the kinetic energy, whereas the local-density profile is 
completely flat.^ 



Figure 3: (Color online) (a) Charge density wave (CDW) and (b) spin-Peierls (SP) phases of SU(4) Hubbard model. In 
both SU(4)-singlet phases, translation symmetry is broken. In CDW [(a)], W-mers {N = 4, here), that are each SU(W)- 
singlet, form (period-2) crystalline structures. In SP [(b)], a six-dimensional representation (6) is formed at each site and 
then it is combined with another 6 on the adjacent site to form an SU(4)-singlet. 


2.3.4. Other commensurate fillings 

The nature of S\J(N) Mott-insulating phases for general commensurate fillings n = p/q 
(p and q being relatively primes) has been investigated with combined use of bosonization and 
DMRG simulations [80]. For N = 2, there is no Mott transition and the charge sector remains 
gapless for any commensurate fillings other than half-filling [41] . 

Again, for N > 2, the physics turns out to be much richer. If q > N, umklapp processes 
are irrelevant and a metallic V-component Luttinger-liquid phase is stabilized, which has N 
gapless degrees of freedom [1 for charge and N — 1 for S\J{N) spin] and hence the central 
charge c = N. When q — N, spin-charge separation occurs and a charge gap opens for finite 
(7 > 0; a gapless Mott-insulating phase emerges with — 1 bosonic modes. The physics is 
then quite similar to the one discussed above (Sec .2.2) for the Sutherland model for the 1/N- 
filling (i.e., one atom per site). Last, when q < N, umklapp processes are strongly relevant 


^ Clearly, this mechanism does not work in the case of = odd where a simple uniform Mott insulator with N/2 
particles at each site is impossible. However, preliminary DMRG simulations showed that we still have an SU(At)-singlet 
dimerized phase with uniform charge distribution. 
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and couple the charge and the spin degrees of freedom. As a consequence, fully gapped Mott- 
insulating phases are formed which spontaneously break the one-site translation symmetry. In 
particular, bond-ordered (dimerized, trimerized, or tetramerized) phases are found depending on 
the filling [80]. For instance, a trimerized phase with a three-fold ground-state degeneracy can 
be stabilized in the SU(6) Hubbard model with two particles per site (n = 1/3). In Fig. 4 we 
provide complementary DMRG data showing the emergence of this phase. 



Figure 4: (Color online) Local density rii and kinetic bond energy for the fermionic SU(6) model with two particles per 
site at f//f = 8 on a L = 72 chain computed by DMRG. There is a strong trimerization pattern with a three-periodicity. 


For filling n = m/N (m = 1,..., iV — 1), a Mott insulator with m atoms per site is formed 
in the large-f/ limit. Then, one can perform a strong-coupling expansion to derive an effective 
spin model for the remaining SU(iV) low-energy degrees of freedom. The resulting magnet takes 
the form of the SU(A^) antiferromagnetic Heisenberg chain with the Hamiltonian [32]: 

H = (24) 

i A=1 

where Sf is the SU(A^) spin operators at site i which transform in the antisymmetric m-tensor 
representation of SU(IV): 

(25) 

For TO = 1 (TV-dimensional fundamental representation) and m = N — 1 (its conjugate), one 
recovers the Sutherland model (19) with a gapless behavior described by SU(A^)i WZW CFT 
[32]. The physical properties of the model (24) for other values of to have been investigated 
by a CFT approach [32, 94] and variational QMC calculations [95]; when to and N have no 
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common divisor, a gapless SU(A^)i WZW quantum criticality emerges.® If N is divisible by 
m (i.e., N = mp), on the other hand, either a fully gapped bond-ordered phase with p-fold 
ground-state degeneracy or the SU(iV)i quantum critical phase realizes depending on whether 
the underlying umklapp operator is relevant or (marginally) irrelevant [32, 95]. The latter case 
corresponds to situations when A^(to—1) > 2to^ [32]. These predictions have been checked with 
a great accuracy recently by variational QMC calculations [95]. When m and N have a common 
divisor, one observes that these results are not in agreement with the SU(iV) generalization of the 
Haldane conjecture proposed in Ref. [96] (see Sec. 3.2.4 for more details of the conjecture). 


filling n 

N = 2 

iV > 3 

incommensurate 

ClSl 

C1S(7V- 1) fort/ > 0, 

CISO for C/ < 0 (MS or CDW) 

1/N 

COSlforC/ > 0, 
CISO for {7 < 0 

C1S(W- 1) forO < C/ < t/c, 

C0S(iV - 1) for U^<U 

1/2 (half-filling) 

same as n = 1 /N 

COSO [dimerized (SP) for U > 0, 

period-2 CDW for U < 0] 



C0S[N — 1) or p-merization 

m/N (m = 1,..., V — 1) 

same as n = 1 /N 

when N = mp, 

C0S(iV — 1) when N and m coprime 



C1S(V — 1) when q > N, 

generic p/q 
{p, q: coprime) 

ClSl 

C0S(iV - 1) when q = N (U > Uc), 
COSO with broken translation 
orC0S(iV- 1) when (7 < N 


Table 1: Phases of single-band SU(A^) Hubbard chain (18) for various fillings. The notation CmcSms (me = 0,1, 
ms = 0,..., — 1) denotes a phase with rric (ms) gapless charge [SU(A^)] degrees of freedom. For incommensurate 

fillings, we have only one gapless charge degree of freedom (CISO) and leading instability is molecular superfluid (MS) 
of A^-mers at low densities or 2fcF-CDW for higher densities (see, e.g.. Fig. 1). The low-energy physics of the phase 
C0${N — 1) is described by the Sutherland model (19) or level-1 SU(A^) WZW CFT. The ‘SP’ and ‘period-2 CDW’ 
states at half-filling are illustrated in Fig. 3. 


3. Two-orbital fermionic Hubbard models 

So far, we have described the phases of the S\J(N) fermions that appear when no additional 
degrees of freedom (e.g., the ground and the excited atomic states of alkaline-earth fermionic 
atoms) are taken into account. However, as we have mentioned in Sec. 1, we can incorporate 
an additional degrees of freedom (“orbital”) into the system by taking into account another set 
of states with SU(iV) symmetry (e.g., the metastable state ‘e’ of alkaline-earth atoms). Due 
to the interplay of S\J(N) and orbital, we have much richer phase diagrams. In this section, we 
consider the one-dimensional SU(iV) fermions having two orbital degrees of freedom and map 
out the ground-state phases. 


®From the effective-field-theory point of view, whether N and m have a common divisor or not affects the selection 
rule for the relevant perturbations allowed at the level-1 S\J{N) WZW fixed point and thereby governs the ground-state 
properties [32]. 
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3.1. Various physical realizations 

In this review, we describe two different ways to introduce additional “orbital” degree of 
freedom. One is to use two atomic states, the ground state ‘g" (^^o) and a metastable excited 
state ‘e’ (^Pq) of the alkaline-earth fermions[17]. Although we use the terminology ‘orbital’ 
here, it is in fact related to an additional internal degree of freedom (g and e) and has nothing 
to do with the real orbital. The highly suppressed internal conversion ^ ^Pq guarantees 
separate conservation of the number of fermions in the g and e states leading to nearly perfect 
U(l) symmetry in the orbital sector. 

The other uses the two degenerate p-bands of a ID optical lattice. Let us consider a ID 
optical lattice (running in the ^-direction) with moderate strength of the harmonic confining 
potential V±(x, y) = -f y^) in the direction (i.e. xy) perpendicular to the chain. The 

two degenerate p-bands are formed by the two degenerate first excited states of the above two- 
dimensional harmonic oscillator and we use them to introduce the orbital degree of freedom. In 
this scheme, the two orbital states originate from the two symmetry-related excited states with 
specific spatial structures. For this reason, we obtain a slightly different effective Hamiltonian for 
the latter case (in particular, we have much less degree of freedom in the effective Hamiltonian). 


3.1.1. g-e model 

Since both g (^S'o) and e (^Pq) states have vanishing total electron angular momentum, the 
same mechanism as in Sec. 2.1.1 leads to the SU(A^)-symmetry in the two-body scattering pro¬ 
cesses (i.e., for all combinations g-g, e-e, and g-e). Even if the scattering lengths for the two 
atoms in a given fixed combination do not depend on the total /, they may be different for dif¬ 
ferent combinations of the atomic states (e.g., g-g and e-e). In general, the scattering length may 
differ for the following four combinations of the two colliding particles: 


Igg) = \g)i\g)2 , |ee) = |e)i|e)2 , 

= ^(lff)i|e)2 ± |e)i|p)2) ■ 

Then, the strength of the interaction is determined by the four s-wave scattering lengths ax 
(X = gg^ ee, ge+, ge~) as 

gx = ~^°-x ■ (27) 

The scattering lengths Ugg, Qee and are for two atoms in the electronic states \gg), |ee) and 
\ge^), respectively. The known values of the scattering length ax for and ®^Sr are summa¬ 
rized in Table 2. In stark contrast to alkali-metal atoms, magnetic Feshbach resonance cannot be 
used, due to the vanishing electron spin, to tune these scattering lengths in alkaline-earth atoms. 
However, quite interestingly, orbital Feshbach resonances have been proposed theoretically [97] 
and observed experimentally in Yb [98, 99], thus allowing to change a^g with a magnetic field. 

When the two-body interaction V is independent of the nuclear spin (this is the case to a 
good approximation in alkaline-earth atoms), the most general form of V may be given by the 
following contact interaction: 


V(r — r') 

= {gggigg}(ggi + gee\ee){ee\+ g+^\ge+){ge+\ + gg^\ge~){ge~\] 5 {r - v') . 
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173Yb (/ = 5/2) 

«7Sr (/ = 9/2) 


10.55 [nm] (Ref. [100]) 

5.09 [nm] (Refs. [101, 102]) 

^ee 

16.2 ± 0.55 [nm] (Ref. [27]) 

9.31 ± 0.58 [nm] (Ref. [26]) 


1 

f 175 ± 16 [nm] (Ref. [28]) 




115±10[nm] (Ref. [27]) 

8.94 ± 0.42 [nm] (Ref. [26]) 


1 

l99±2[nm] (Ref. [98]) 


age 

11.6 ±0.1 [nm] (Refs. [27,98]) 

3.60 ± 1.16 [nm] (Ref. [26]) 

^qe ^qe 

111 [nm] (Ref. [99]) 



Table 2: Values of scattering length ax = QQ, ee, ge^, ge ) known from experiments. 


Note that all the four couplings Qgg, gee, and are independent of the nuclear-spin states of the 
colliding atoms (see the discussion of Sec. 2.1.1). From (28), the orbital-dependent part of the 
matrix elements is calculated easily: 


(u, 6|f^|77r, 7l) — ggg^a,g^b,g^m,g^n,g gee^a,e^b,e^m^e^n, 

1 


2 ^9^ i^a,g^b,e^m,g^n,e “t” ^a^e^b,g^m,e^n,g) 

“f “ ^9ge 9ge') i^a,g^b,e^m,e^n,g ^ ^(t* t* ) 


(29) 


{a,b,m,n = e,g) . 

The derivation of the lattice Hamiltonian for the case of two orbitals closely follows that 
described in Sec. 2.1.2 for the single-band case except that now we have two species of fermions 
and use the corresponding Wannier functions: 




R 


R nGbands 


(n) ^ , 1 

C'r 




VN, 


cell 


E 


R n 

{m = e,g] a=l,...,N) . 


(30) 


We plug these operators to rewrite the two-body interaction 

1 , f f 

2 E E J dz J dz'cle,{z)clp{z'){a,b\v\m,n)cn/ 3 {z')crnaiz) . (31) 

06,/3—1 a,b,m^n—e,g 

In contrast to the single-band case where we have obtained only the Hubbard-(7 interactions, we 
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now have the following four different types of interactions: 


1 

2 


E 


E 


Tr /ni,n2;n3,"4 \ ("l)t (" 2 )t (" 3 ) (" 4 ) 

vHi,i?25-R4 9<^-,Rig&,R2 g^iRz 9 cic,Ra 


1 ^ 

■ 2 E E 

q,/3 —1 {Ri,ni} 


ni,n2;n3,n4 

Ri ,R2‘,R3yR4 


)c 


ea,Ri 


„("2)t „(n3) (" 4 ) 

^e^,i ?2 ^ea,i?4 


1 _ 

_ Tr+/ni,n2;n3,n4 \ „(’4l)t „(’42)t „("3) „("4) 

2 / ^ ^ge\Ri,R2-,R3,Ri) ^ga,Riefi,R2^eP,R3^ga,Ri 


q,/ 3—1 {it^,ni} 


+ 


1 

2 




Y r„( 


Rl ,i?2 ;i?3 ,i?4 


'1 „('*l)t „(n2)t J’T’3) (" 4 ) 

/ ^ga,Ri^el3,R2^9^,R3^ecx,R4 ‘ 


{Ri,ni} 


In the above, the interactions are given in terms of the Wannier basis as 


(32) 


T/ ;'^3 5'^4 

mm \Ri ,i?2 ;i?3 ,i?4 


) = 9aa f dz {z)wi'^R, (z) (m = g, e) 


J/+ {ni,n2-,n3,n4 

^ge\Ri ,R2-,R3,R4 


) = i9ge + 9. 






t/“ ('"1’”2;"3,™4 \ _ ( + 

''ge\R^,R2-,R3,Ri' — ^dge 



dz W 


(ni)* 

g,Ri 


{z)w 


{ 712 )* 

e,R2 


(z)w 


("3) 

g,^! 


(z)w 


(ni) 

e,i?4 


( 2 ), 


(33) 


and the coupling constants ggg, g^e, and g'^^ are calculated from the scattering lengths as Eq. (27). 
The first two terms in Eq. (32) describe the density-density interactions between the fermions 
in the same orbital, while the third and the fourth ones correspond to fermions from different 
orbitals. Specifically, the third one is just the density-density interaction of a pair of fermions 
on different orbitals and the last one is the orbital exchange interaction (or, the Hund coupling). 
Here it should be noted that there is no special relation among the four couplings Vgg, I4e, 

as the Wannier functions for the two orbitals ^ (,z) are not symmetry-related to each other^. 
As before, we restrict ourselves only to the same band ni = 712 = ns = 714 = uq and keep only 
the onsite terms Ri = R 2 = R 3 = R 4 = i to obtain the following Hamiltonian (g-e model) [17]: 


N 

T~Lg-e=— (^^Ina. i^ma, i+1 + H.C.) 

m—g,e i a — 1 

- E - 1) (34) 

m—g,e i m—g,e i 

+ v E Rg,i'Re,i~\~ i^e/3, 2^3/5^ 

i i,Qt0 


where the index a labels the nuclear-spin multiplet and the orbital indices m = g and e label the 
two atomic states ^Sq and respectively. We have also introduced the number of the fermion 


^This is not the case for the p-band model where the two orbitals are related to each other by C 4 -symmetry. 
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m at each site 


N 

‘^m,i = ^ ^ — 5 ;^) ( 35 ) 

a=l 


and suppressed the common band index no- The coupling constants are given in terms of Vgg, 
Vee, and as [see Eq. (33)] 


Urr 




T/ _ {no,nQ-,nQ,no \ j/Q-e _ xr—/no,no;no,no \ 

^ — ^g€\Ri,Ri;Ri,Ri)^ ^ex “ ^ qe\R^ .Ri :R^ .RJ 


(36a) 


while the hopping tm and the chemical potential are given by the Bloch energy as 


tmU-f) = (m = g,e) . 

-‘''cell , 


(36b) 


In order to understand the physical processes contained in this Hamiltonian, it is helpful to 
represent it as two coupled (single-band) SU(A^) Hubbard chains (see Fig. 5). On each chain, 
we have the standard hopping tm along each chain and the Hubbard-type interaction Umm, and 
the two are coupled to each other by the nearest-neighbor Coulomb interaction V and the g-e 
exchange process 14® The ‘hopping’ between the chains does not exist as the transition g e 
is strongly suppressed. Therefore, on top of the obvious SU(A^)-symmetry, the Hamiltonian is 
invariant under the following U(l)-symmetry 


i ' t e i ; 


V 


(37) 


This is a consequence of the fact that the total fermion numbers for g and e are conserved sepa¬ 
rately.^ 



Figure 5: (Color online) Two-leg ladder representation of the g-e Hamiltonian (34). 


To understand the global phase structure, it is useful to rewrite the exchange interaction V^x ^ 
in two different ways. First, we introduce the second-quantized SU(A^) generators of each orbital 
as 

N 

Sm,i= cY,i{S^)apCm 0 ,i {m = g,e). (38) 

a,/3—1 


“This breaks down when there is transition (i.e., ‘hopping’) between the two atomic states g and e. 
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If we normalize the SU(A^) generators {5'"^} as® 


Tt{S^S^) = 5^^ , 


(39) 


they satisfy the following identity: 

^ (iS"^)a/3(‘5"^)7<5 = f — ^<5a/3^75 j • (40) 

A=1 ^ ' 

Then, it is straightforward to show that the orbital-exchange interaction 14|'® can be written as 
the ‘Hund coupling’ between the SU(iV) ‘spins’ of the two orbitals 

\ 1 

^^ I ^ ^g,i^e,i I ~ J^'^ng^ine,i ■ (41) 

i,a0 i \ A=1 J i 

The fermionic anti-commutation is crucial in obtaining the minus sign in front of the Hund 
coupling. The above expression enables us to rewrite the original g-e Hamiltonian (34) as [103] 


n 


g-e 


N 

“ X! X! X! ( •'ma, i-l-1 + H.C.^ 

i m—g,e oc—l 

m—g,e i i m—g,e 





Hund 


(42) 


From this, one readily sees that positive 14® ® leads to ferromagnetic coupling between the two 
SU(A^) spins on the g and e orbitals thereby maximizing the SU(iV) ‘spin’ at each site. 

To derive the third form, we first introduce the orbital pseudo-spin 


N N 

2 ’.“ = _ (f CF°- c ■ = T® 

q;=1 m,n—e,g a=l 

(a = X, y, z; m,n = g, e; (t“ :Pauli matrices) . 


(43) 


Clearly, the U(l)-symmetry (37) is generated by Tf and, in the following, we call it U(l)o. 


^This corresponds to, e.g., using the SU(2) generators (t“/\/ 2 (a = x, y, z) instead of the standard ones t7“/2. 
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Using these pseudo-spin operators, we obtain the third form of the g-e Hamiltonian [103] 


N 

T-Lg-e= — '^2i ('^Ina. jCmg, i +1 + H.C.) 

i m—g,e a—1 

- E (M^'"^+3Kr/4)E^™.^ 

m—g,e i 

Umm — Vji‘^/2 

E / ^ ^ ^ 1) 

i m—g^e 

+ {V + Vr/2) rig, .n,,. + E(^*)' ' 


Hund 


(44) 


This expression, which is equivalent to (42), represents the exchange interaction 141'^ in terms 
of the orbital SU(2) pseudo-spin T. It is important to note that the sign of 14l'® is opposite to 
the one in Eq. (42); positive (negative) 14® tends to quench (maximize) the orbital pseudo-spin 
and maximize (quench) the SU(A^) spin. This dual nature of the orbital pseudo-spin T and the 
SU(A^) spin is the key to understand the structure of the phase diagram and we will come back 
to this point in Sec. 3.2.1. 



Figure 6: (Color online) Typical strong-coupling ground states for SU(4) g-e Hamiltonian (42) with tm = 0. The SU(A^) 
‘magnetic moment’ that appears at each site in the Mott-insulating limit is shown by Young diagrams. For ^ < 0, 
SU(A^)-spin at each site is quenched [see Eq. (42)], while orbital pseudo-spin T is quenched (and hence SU(A^) spin is 
maximized) for 14 ? ^ > 0. 


Before concluding this subsection, let us give another useful form of the g-e Hamiltonian 
[103]: 


N 

^5-e = ~ E/ E/ E! (rilna. jl^rna, i+1 + H.C.) 

i m—g,eoL—l 

~ 2 4" fTg) E! rii + E! rii 

i i 

+ ^E {(E)" + (E)"} + - {pg - Me)EE 

i i i 

+ t^diff E ) 


(45) 
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where a set of the new coupling constants are given in terms of those in Eq. (34) 

U = \{Uag + U,e + 2V), t/diff = ^{Ugg - Uee), 

J=vr, J. = ^{Ugg + Uee-2V), (46) 

Ms = + Ugg + vr). Me = + C/ee + KD ■ 

The model (45) with Ugg = Uee, M^®^ = M^*^^ dubbed the generalized Hund model and has 
been studied extensively for TV = 2 in the cold-fermion context [104, 105]. It is obvious that, 
when tg = tf,, J = Jz = Udiff = 0, p-g = pe. the Hamiltonian Hg.e is U(2TV)-invariant and the 
orbital part (J and J^) breaks it down to the generic symmetry U(l)(,xSU(TV)j,xU(l)o: 


U(2iV) U(l)eXSU(iV),xSU(2)o 

C^diff=0 

-^^^U(l),xSU(TV)^xU(l)„ . 

or (7diff#0 


(47) 


In Sec. 3.4, we use the form (45) of the g-e model and the parametrization tg = J = Jz = 
t^diff = 0, fig = He to map out its phases. 


3.1.2. p-band model 

In the second scheme, we use only the ^-state (^^i) and, to implement the orbital degree 
of freedom, introduce the two degenerate p-bands of a ID optical lattice. Let us consider a 
ID optical lattice (running in the ^-direction) with moderate strength of (harmonic) confining 
potential V±{x,y) = ^Muj^y{x^ + y^) in the direction (i.e. xy) perpendicular to the chain. 
Then, the single-particle part of the Hamiltonian reads as 

"" = (48) 

= 'H//(z) +n_L{x,y) , 


where Iper(^:) is a periodic potential that introduces a lattice structure along the chain (i.e. z) di¬ 
rection. If the chain is infinite in the z-direction, the single-particle state is given by the following 
Bloch function: 


Ja:, p, z) = , 


(49) 


where the two functions (z) and (j)n^ (a:, y) respectively are the eigenfunctions of'H//{z) 

and j_ {x, y) (see Fig. 7): 


n//{z)vtSz) = e<^'^Hkz)pt^{z) (50a) 

T~Lx{x,y')(j>{n^^ny){x,y) = e{n^,ny)4'(n^,ny)(.X,y) . (50b) 


with 


^(n^.riy) — (^(c 4“ T 1) kWxy {j^xi tty — O 5 1; 2, . . .) . (50c) 

This implies that due to the motion perpendicular to the chain, each one-dimensional Bloch band 
specified by n splits into subbands labeled by {ny;,ny): 
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(51) 




We call, for a given main band index n, the subbands with (n^;, Uy) = (0, 0), (1, 0), and (0,1) 
as ‘s’, ‘px’ and ‘Py \ respectively. 



-3 -2 -10 1 2 3 -3 -2 -1012 3 -3 -2 -10 1 2 3 


x/xo 


x/xo 


x/xo 


Figure 7: (Color online) Contour plots of squared wave functions \ (l)nj;,ny P for three orbitals {rix^riy) = (0,0), (1,0) 
and (0,1). From Ref. [103] 


Now let us suppose that only the n = 0 bands are occupied, and that, among them, the lowest 
one (the s-band) is completely filled. Then, it is legitimate to keep the next two bands and 
Py in the effective Hamiltonian in describing the low-energy physics [106, 107]. The argument 
goes almost in the same way except that now we have to take into account the motion in the 
perpendicular (i.e., xy) directions. Now the Wannier function 


w 


(n) 
(n^,ny);R'~ 


(x,y,z) = —)(a;,y)^ 

V^cell 




(z) 


(52) 


(H labels the center of the Wannier function and iVceii is the number of unit cells) is used instead 
to expand the Fermi operators; 




I ,a,Ii 


(a=l,...,N). (53) 


R n=bands (rix ,ny) 


We plug them into the two-body (contact) interaction for atoms in the g state 

\ X! fdr fdr'cl{r)cl{r')V{r-r')ci3(r')cair) (54) 

q./3=1 


with 

V[v-v') = ggg5{v-v') . (55) 

Retaining only the terms with Ri = R 2 = R 3 = R 4 = i, ni = n 2 = = 0, 

{uj;, Uy) = (1, 0), (0,1), we obtain 


N 


2 Uaaaa ,RCa,a,R 

R ct,/3=l Ka=p^,py 

+ Uabba <^a,a,R^b,P,R^b,P,RCa,a,R 


a^b 
-Px ,Py 


a^b 

^Px,Py 


(56) 


^abab cl^^^^iicl^^ j^Ca,p,RCb,a,R\ , 
a^b ) 
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where the superscript ‘(0)’ for the fermion operators of the lowest Bloch band has been sup¬ 
pressed. We have also introduced a short-hand notation m = Px,Py with p^ = {ux, Uy) = (1,0) 
and py = {ux^riy) = ( 0 , 1 ). 

As the Wannier functions are related to each other by the C 4 -symmetry, the coupling 
constants Uabcd defined by 

Uabcd = Pgg {a,b,c,d = Px,Py) (57) 

obey the following relation: 

“ ^PyPvPvPv = 

^PmPmPyPy ~ ^PyPyPxPx ~ ^PxPyPyP:^ ~ ^PjjPxPxPh = ^PmPyP^Py ~ ^PyPxPyPx^ ( 5 ^) 

= U2 . 

Therefore, in contrast to the case of the g-e model (34), we have only two independent couplings 
Ui and U 2 for any C' 4 -symmetric potentials. In fact, further simplification occurs for axially- 
symmetric potentials like the harmonic one used here. In these cases, {x, y) in Eq. (52) is 

replaced by 

f{r) cos 9 for px , f{r) sin 9 for py [r = \/x'^ + y'^) (59) 

with some real-valued function /(r) depending the actual choice of the potential. Carrying out 
the integration in Eq. (57), we obtain 


Ui = ^Ir,x, U2='^Ir,z, (60) 

where Xr^z is the result of the integration over r and z that is common to all Uabcd- This implies 
that {Ui, U 2 ) are bound to satisfy the relation 


Ui = W 2 (61) 

for any axially-symmetric potentials V± (x, y) = fi\/x‘^ -I- j/^). That is, as far as we use axially- 
symmetric traps, there is only one free parameter for the interactions. One possible way to 
deviate from the line Ui = 3 U 2 is to use such (C' 4 -symmetric) anharmonic potentials as [103]: 

V±{x,y) = ^Mujly{x‘^+ y^) + ^f3{x'^+ y‘^) (/3 > 0) (62) 

(see Eig. 3 of Ref. [103] for the ratio U 1 /U 2 obtained for the above El). 

As the single-particle energy ^ )(^z) [Eq. (51)] is the same for the two p-bands, we 
have the same hopping amplitude and the chemical potential for px and py. Summing up all 
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these, we obtain the following form of the p-band model [103, 106, 107]: 


^p-band — ip ^ ^ ~i~ H-C-) 

i m=p^ ,py 

Mp 'y ^ y ^ “t” 2^^ y ' ^ ^ 2(^771, 2 1) 

2 m=p^,py i m=p^,py 


(63) 


E 


^lna,i^lnp,i^ri,p,iCna,i , 


2 q:,/3=1 m^n 

—Px,Py 


where the hopping tp and the chemical potential pp are given by the single-particle energy 

=- ^ E ^ =- 7^ E ^fri) (y) ■ ( 64 ) 


cell 


cell 


The last term of Eq. (63) comes from the pair-hopping between the two orbitals, which is not 
allowed for the setting of the g-e model, and breaks U(l)o-symmetry in general, while the other 
five terms already existed in the g-e model (34). In fact, except for the last term, "Hp-band coincides 
with the Hamiltonian Hg.e [(34)] after the following identification (see Fig. 8) 


tg — tg — tp, — y — i^p ^ 

Ugg =Uee = Ul, V = U2, Kt" = U 2 ■ 


(65) 


Due to the pair-hopping between the two orbitals, the p-band Hamiltonian (63) appears to break 
the U(l)o-symmetry. To see this, we rewrite (63) using the orbital pseudo-spin: 


n 


p-band 


= — t 


E E ( ^ma, i^ma, 2+1 + H.C.^ 

2 irL—P:c,Py 

- mpE^® 4- + U 2 ) y 

i 2 

+ ^{2(72(T)+ + ((/i-C/2 )(T’+} , 


( 66 ) 


which is to be compared with Eq. (45). As can be easily seen, the last two terms break U(l)o in 
general. Therefore, the generic symmetry of the p-band model is 


U(l)eXSU(Af),xZ 2 ,o . (67) 

However, this is not always the case and, in fact, there is a hidden U(l)o symmetry in the case of 
axially-symmetric traps. In fact, if one plugs into (66) the relation Ui = ZU 2 that holds for any 
axially-symmetric V 2 _{x, y) [see Eq. (61)], the orbital part assumes a fully U(l)-symmetric form: 
2 U 2 -I- ( 17 )^} and Hp-band reduces to a special case of Hg.e [Eq. (45)] (with pg = pe, 

Udis = 0) after the redefinition: -O- Tf, T/ —)■ —Tf}^ 


'^This is in a sense an artifact of the choice of the basis (p^ and py). In fact, if we had chosen the orbital-angular- 
momentum (along the 2 -axis) basis, the U(l)o-symmetry would have been explicit. 


25 





A few remarks are in order here about the special points where great simplification or en¬ 
hanced symmetries emerge. When U 2 = 0, the p-band model decouples into two non-interacting 
copies of SU(A^) Hubbard chains and we can borrow the results from Sec. 2. Moreover, along the 
line Ui = U 2 , the p-band model (66) reduces to the above U 2 = 0 case (two decoupled chains) 
after the redefinition o Tf, Ui —> Ui/2. Finally, in the N = 2 case, the p-band model can 
be recast in the following form [107]: 

'Hp-hand = — tp + l +H.C.) 

'm=p^ ,py j 

2 / \ ^ 

~ ~ ( 68 ) 

i i rn—px ,Py 

- 2 U 2 E 2t72 YyKp^,,-kp^,, 


using the following two commuting sets of spin and pseudo-spin operators {.§“} and {7^“} 
[108]: 


{a = x,y,z; m=p^,Py) 
km,i — (~1) -^m,i = (~1) 

km,i = 2 ~ 1) = ~ 1) • 


(69a) 

(69b) 


Note that the set of operators K generates the so-called 77-SU(2)-symmetry in the charge sec¬ 
tor. From this, we can see that the N = 2 p-band model at half-filling enjoys an enlarged 
SU(2)5,pi„ X SU(2)(,[jafge ~ SO(4) symmetry for all Ui, U 2 which stems from an additional SU(2) 
symmetry for the charge degrees of freedom at half-filling [107]. For more details about the 
SO(4) symmetry in the Hubbard model, see section 2.2 of Ref [60]. 



Px"00 




U 2 


Figure 8: (Color online) The two-leg ladder representation of the p-band model (63). On top of the interactions included 
already in the g-e model, pair-hopping processes between the two orbitals are allowed. 


3.2. Results known for Mott-insulating limits 

As in the usual Hubbard Hamiltonian, when we have an integer number of particles at each 
site, the system becomes a Mott insulator in the limit 

U =^{Ugg+U^, + 2V)^<X) (70) 

[see Eq. (45)]. Depending on N and the number of fermions at each site, we have a variety of 
SU(A^) spins emerging at each site in the Mott-insulating limits (see, e.g.. Fig. 6). 
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3.2.1. Low-energy degrees of freedom in Mott phases 

To find the SU(-/V) and the orbital SU(2) content, we first note that when J = VS.'^ = 0, 
Jz = 0, Udiff = 0, the system attains the maximal symmetry U(2A^) [see (47)]. In this limit, 
for a given number of fermions per site (say, n), we can uniquely assign the following U(2A^) 
irreducible representation: 

(0 < n < 2iV) . (71) 

Formally, the ^]^N-ny. states in the above representations may be decomposed in terms of 
the direct product of SU(iV) and the orbital SU(2). This problem is explicitly worked out in 
Ref. [109] (see also Appendix A of Ref. [103] for a concise explanation in the context of the 
SU(A^) Hubbard model). As a simplest example, let us consider the n = 1 case. Then the 2N 
singly-occupied states at each site are decomposed as 

^ = T = l/2). (72) 

SU(2N) SV(N) SU(2) 



Therefore, the low-energy physics of the highly degenerate ground-state manifold is described 
by the composite of the SU(A^) ‘spins’ in the fundamental representation (□) and the orbital 
pseudo-spins T = 1/2 [17]. 

Different results are obtained when we consider the case of half filling. For instance, the 70 
half-filled states of the SU(4) g-e model (i.e., = 4, n = 4) may be decomposed as 



= (20, T = 0) © (15, T = 1) © (1, T = 2) 

(with • being the singlet). Note that the Fermi statistics dictates the possible combinations of 
SU(A^) and orbital. The Hund coupling 

/ n ^-1 \ 

- E E + (density-density interaction) (74) 

i y A=1 J i 


[see Eqs. (42) and (44)] then selects one of the irreducible representations appearing on the right- 
hand side: 



for > 0 
for < 0 . 


(75) 


That is, when the SU(A^) spin is maximized (quenched), the orbital pseudo-spin T is quenched 
(maximized) as is expected from the discussion in Sec. 3.1.1. This dual nature of the SU(A^) and 
the orbital is common to all SU(A^) g-e models with even N. 
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For odd-iV, on the other hand, we have different decompositions. For instance, for TV = 3 
and n = 3 (half-filling), the decomposition in terms of SU(3)xSU(2) reads as 


n = 3 



\ 


SU(2) 


(• > 


\ SU(3) 

(8,r = i/ 2 ) 0 (i,r = 3/2) . 


(76) 


Therefore, the situation on the 14® ® > 0 side is rather different for iV-even and iV-odd. Specifi¬ 
cally, the low-energy degrees of freedom in the half-filled Mott insulator for A^-even are the pure 
SU(A^) spins, while those for A^-odd are the composite of SU(A^) spins [with {N + l)/2 boxes 
in the first column and {N — l)/2 boxes in the second] and orbital pseudo-spins T = 1/2. 


3.2.2. Strong-coupling limits 

In the previous section, we have seen that different low-energy degrees of freedom emerge 
depending on N, filling n(G Z), and the sign of Therefore, when the hopping amplitudes 
vanish (atomic limit), i.e. tg = te = 0 or tp = 0 depending on the model, the ground state 
has extensive degeneracy. This degeneracy is lifted when the hopping tm is taken into account. 
Clearly, changes the fermion number n at each site and non-trivial matrix elements require at 
least processes. 

For the clarity of the argument, we restrict ourselves to the case of half-filling (i.e., the number 
of fermions per site n = N) in the following. Let us begin with the case where N = even, 
Vl"" > 0. Then, in the atomic-limit ground state, we have 


N/2 



(77) 


with the orbital pseudo-spin quenched (T = 0). The resulting effective Hamiltonian is obtained 
by the second-order perturbation and reads as follows [103] 


N^-l 

7fsu(7V) = <7su(Ar) ^ Sf" 5/^2 + const. , 

A=1 


where the exchange coupling Jsv(N) is A^-independent; 


J; 


ti 


ti 


SV(N) = 


U -f Udiff + J + -0 U — [/diff + J + -0 


(78) 


(79) 


and the SU(AT) spins {5^} transform under the representation (77) [see Eq. (75)]. As T is 
quenched, we obtain a pure ‘spin’ Hamiltonian (as we have seen above, this is not the case 
when N = odd). 

In the case of Tfp-band, the condition I4x > 0 translates to U 2 > 0. However, is not 
conserved in general and we cannot use the argument in Sec. 3.2.1 as it is. However, we found 
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that when Ui > U 2 {> 0) the lowest-energy state is T = 0 singlet that enables us to follow 
exactly the same steps and obtain [103] 


■Hsuttv) = .. ■ (^0) 

In Sec. 3.3, we will show that the SU(iV) spin Hamiltonians (78) and (80) host a topological 
phase in quite a large region of the parameter space. 

When 14x ® < 0 (or U 2 < 0), on the other hand, the low-energy degree of freedom is the 
orbital pseudo-spin transforming like a T = N/2 spin [Eq. (75)] 



N boxes (T=N/2) 


(81) 


and the SU(7V) ‘spins’ are quenched. The resulting effective Hamiltonian is rather different 
for the g-e Hamiltonian and the p-band Hamiltonian. For the g-e Hamiltonian, we obtain the 
following T = N/2 Hamiltonian; 


{j.y - (J - J.)iTtr} 

i 

i 


(82) 


with the easy-axis exchange coupling [103] 

^tgte _ 

U+\J\ (iV+i)} 

2U^+t^] 

n{u + \j\{n + ^)} 

In the ideal situation where the two orbitals are symmetric, the model reduces to the spin T = 
N/2 Heisenberg model ( Jxy = Jz) with the single-ion D-term, for which much is known (see, 
e.g. Refs. [110-112] and references cited therein). When \ Jz\ <C | J| (J < 0), the system has 
easy-plane anisotropy and, if N is even^' and | J| is large enough compared with J^y and Jz, the 
ground state is in the so-called large-Z? phase, which, in its extreme case, is a product of T* = 0 
states. In the fermionic language, it is a state where there are N/2 fermions on each orbital (and 
hence T/ = 0) and these N fermions form total SU(W)-singlet at each site. 

When Jz takes a large negative value (this is the case typically for large V), the ZJ-term is 
of easy-axis-type and the ground state is Ising-like: = • • • , —N/2, N/2, —N/2, N/2, ■ ■ ■. 

This is the 2fcF(= tt) density wave of the orbital pseudo-spin and we call it orbital density wave 
(ODW). See Sec. 3.4 for the location of the phase. 

In the case of the p-band model, physics of the orbital sector is even more involved. Since the 
condition 14f’® = J < 0 translates to f /2 < 0 in the p-band model [see Eq. (65)], the interaction 



* ’ When N is odd, the large-D limit is not very trivial and the effective Hamiltonian is given by the S = 1/2 XXZ 
chain. 
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is attractive J7i + C /2 < 0 in the physical region C7i ~ 3 C /2 and we have to take into account 
several different values of rii. Therefore, the physics here is quite different from the usual Mott 
physics described above. For instance, at p = —N\Ui + C/ 2 I, we have two degenerate S\J(N)- 
singlet states = 0 (T = 0) and Ui = 2N (T = 0) which feel a repulsive interaction coming 
from f^-processes thereby stabilizing the 2fcF-CDW*^ in a region around the line Ui = 3C72(< 0) 
for iV > 3 (see, e.g.. Fig. 22 of Ref. [103]). 

Last, the case with = J = 0 is slightly different. This region is most conveniently 
investigated in the limit U = Ugg = Uee = V ^ t, Then the g-e model (45) 

reduces to the single-band SU(2A^) Hubbard model at half-filling discussed in Sec. 2.3.3; the 
ground state is a dimerized SP state (see Fig. 3). This highlights the importance of the orbital 
degree of freedom and the exchange interaction between the two orbital states in realizing 
topological phases. 

3.2.3. Sigma-model mapping 

Some insights into the nature of the ground states may be gained by mapping the problems 
onto the sigma model with the 0-term. Generalizing the semi-classical sigma-model mapping 
[57, 1 13] a la Haldane of the usual spin chains. Read and Sachdev [114-116] considered a family 
of SU(iV) ‘spin systems’ where an irreducible representation TZ (specified by a rectangular Young 
diagram with m rows and columns) and its conjugate TZ are assigned respectively on the two 
sublattices of a //-dimensional bipartite lattice. 


7 ^ = 



ric 



(84) 


If we represent the state (84) in terms of tlc copies of SU(/V) fermions, these states are in fact 
invariant under \J(N — m)x\J(m). Therefore, the coherent state, generated from the reference 
state (84) by applying any elements of U(A^), represents ;j(jv u(m) "This is closely 

parallel to that the Bloch spin coherent state generated from the spin state polarized in the z- 
direction is isomorphic to ~ u(px^u(i) “ ~ expectation value {S) 

of spins in the usual spin path integral is now played by an A^x hermitian matrix Q {Q^ = 1). 
Since the number of columns ric controls the semi-classical limit as the spin S does in the SU(2) 
case, we may expect that an S\J(N) analogue of Neel ordering occurs in the limit ric — 00 . The 
order parameter is the staggered component H of the matrix Q and its low-energy fluctuations 
are governed by the following (hnite-temperature) effective action [115]: 




rTr{a^H(T, r)}^ 



(85) 


where the two phenomenological parameters Xs and Ps{v)fj,jy respectively are the transverse sus- 


^^The basic mechanism underlying the CDW here is the same as the CDW in the case of the half-filled single-band 
SU(A^) Hubbard chain discussed in Sec. 2.3.3. 
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ceptibility and the spin stiffness'^ [115, 117]. In one dimension (D = 1), this action is supple¬ 
mented by the 6*-term'"^ associated with 112 (^ 'u(n -’^■^xu(m) ) ~ 

0 rPH r 

= dr J dxe^^TT{n{T,x)df,n{T,x)dM{T,x)} 

— “bl^ €)top — TIq'k') . 

Following the same line of argument as the one first used by Haldane [118, 119], we may con¬ 
clude, for D — 1, that when the number of columns ric is even for which the 0-term is irrelevant 
to bulk properties, the system is in a featureless gapped phase [115, 116]. On the other hand, 
when Tic is odd, the system spontaneously dimerizes. In the context of our SU(iV) Hubbard 
model, the case ric = 1 corresponds to the Mott-insulating phase discussed in Sec. 2.3.3 (the 
single-band SU(A^) Hubbard model at half-filling)'^ and the model with ric = 2 describes the 
half-filled 2-orbital system (Sec. 3) deep in the Mott-insulating region. For instance, the appear¬ 
ance of the SP phase in the single-band model [45, 82] discussed in Sec. 2.3.3 is consistent with 
the above field-theory prediction. In Sec. 3.3, we will argue that the featureless phase predicted 
by the model (85) is in fact topological by explicitly constructing the model wave function. 

3.2.4. Lieb-Schultz-Mattis argument and generalized Haldane conjecture 

The simple argument first used by Lieb, Schultz, and Mattis [120] is quite general but some¬ 
times able to give strong constraints on the nature of the ground state and the low-energy spec¬ 
trum over it. In fact, Affleck and Lieb [121] extended the original argument of Ref. [120] to in¬ 
clude the case of self-conjugate representations of SU(iV). Specifically, they showed for SU(A^) 
{N = even) chains based on the self-conjugate representation with N/2 rows and columns 
that (i) if ric = 1 [anti-symmetric (A^/2)-tensor representation], the finite-size ground state is 
unique. Furthermore, assuming the uniqueness of the ground state, they showed that (ii) for any 
translationally invariant choice of representations'® (i.e., the same representation TZ is assigned 
for all sites), the SU(iV) chain harbors low-lying excitations whose gap is bounded by 1 jL {L be¬ 
ing the system size) provided that the number of boxes riy in the Young diagram representing TZ 
is not divisible by N. In other words, except for the cases of ny = 0 (mod N) [including the one 
with N/2 (N = even) rows and two columns which is relevant to our effective spin model (78)], 
this statement excludes the possibility of gapped topological ground states. Remarkably, this is 
perfectly consistent with the recent group-cohomology classification of the gapped topological 
phases in 1D[122] (see Sec. 3.3.2 for the detail). 

In the context of the SU(A^) (N = even) spin chains based on the representation with N/2 
rows and ric columns (self-conjugate representation; the model (78) obtained in Sec. 3.2.2 corre¬ 
sponds to ric = 2), this means a unique gapless ground state or gapped degenerate ground states 


•^For the nearest-neighbor Heisenberg SU(A^) magnets on a hypercubic lattice in D dimensions, 


I . ^ 2-Dr 

Ps(r)/.iz = —a^ 5 ^^, x^ = 


16D Ja^ 


( 86 ) 


(with ao being the lattice constant). 

^^In contrast to the corresponding expression for the 0(3) non-linear sigma model, an extra factor i is not necessary 
here as the integral itself is pure imaginary. In the case of the 0(3) non-linear sigma model, we express the matrix field 
by a unit vector n as O = n-cr to recover the factor i in the 0-term. 

^^As we are dealing in this paper with the translationally invariant systems, 7^ = 1^, i.e., m = N/2. 

^^The proof of the existence of low-lying states works regardless of whether the ground state is unique or not. However, 
unless the (finite-size) ground state is unique, the proof does not tell anything about the excited states. 
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unless Nxn^./{2N) = 0 (mod 1), i.e., n^. = even. This is consistent with the argument in the 
previous section since the sigma-model mapping tells us that when = UcTt = 0 (mod 2tt) 
the ground state is unique and gapped.'^ 

There is an SU(4) spin model with two exact charge-conjugation (or translation) breaking 
ground states [123]. According to the Lieb-Schultz-Mattis argument [121], the physical spin 
(ny 7 ^ 0 mod 4) at each site implies that the (thermodynamic-limit) ground state is either gapless 
or (at least two-fold) degenerate. The existence of the two degenerate ground states means that 
the second possibility realizes here. 

All the above statements concern ID systems. In the SU(2) cases, the higher-dimensional 
extension of the Lieb-Schultz-Mattis argument is known [124, 125]. A similar extension for the 
SU(A^) cases would be intriguing in view of the proposal of exotic spin liquids in two dimensions 
[126, 127]. 

Last, we mention a generalization of the Haldane’s conjecture for SU(2) spin chains to SU(A^) 
proposed in Ref. [96] (see also Ref. [128]). Three different cases (type-I, II, and III) have been 
introduced there depending on N and ny- Type I concerns the case when N and ny have no 
common divisor and an SU(A^)i quantum criticality with central charge c = N — lis expected. 
When ny is divisible by N (type-II), a Haldane-gap phase is expected in general (note that this 
case is not covered by the theorem of Ref [121]). The last category (type-III) corresponds to 
the case when ny and N have a common divisor different from N. For short-range interactions, 
an SU(W)i quantum critical behavior is expected. In this respect, some DMRG calculations for 
A^ = 4 and ny = 2, i.e. type III behavior, have found an SU(4)i criticality [96]. As already 
noticed in Sec. 2.3.4, however, this conjecture is not consistent with the known results that SU(A^) 
Heisenberg spin chain in self-conjugate antisymmetric representations (i.e., type III case) have a 
spectral and a ground-state degeneracy. 

3.3. Symmetry-protected topological phases 

In this section, we try to characterize the nature of the ground state of the S\J{N) spin chain 
(78) obtained deep in the Mott-insulating phase. Specifically, we show that the ground state of 
the model (78) shares essentially the same properties with that of the solvable valence bond solids 
(VBS) models in Sec. 3.3.3 and that it is in fact in one of the symmetry-protected topological 
(SPT) phases. The concept of SPT phases is a generalization of topological insulators [129] to 
a class of interacting states of matter that are characterized by short-range entanglement and are 
well-defined only in the presence of certain symmetries (called protecting symmetries) [37, 39]. 
Being topological, this class of topological phases defies the traditional characterization with 
broken symmetry and the corresponding local order parameters; two phases are distinguished 
not by the symmetry they possess but by the presence of quantum phase transitions that separate 
them [39]. One way of distinguishing between topological phases from trivial ones and labeling 
the former is to use the physical edge states. However, this approach is not quite satisfactory 
in the following respects. First, even topologically trivial states may have certain structures 
around the edges of the system, as, e.g., in the spin-2 Heisenberg chain [130, 131]. Second, in 
order to see the edge excitations, it is necessary to consider the excitation spectrum, while the 
topological properties are intrinsic to the ground state itself and should be seen only by examining 
the ground-state wave function. 


*^When ric = odd, on the other hand, the sigma-model argument generically predicts dimerized ground states except 
for N = 2, where we expect the gapless SU(2)i WZW criticality. 
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Recently, the use of the entanglement spectrum, which is defined as the logarithm of the 
spectrum of the reduced density matrix, in characterizing topological phases has been suggested 
in Ref. [132]. This is based on the observation that the entanglement spectrum resembles the 
spectrum of the physical edge excitations. The idea has been successfully applied to various 
ID systems [133-139] and enabled us to characterize topological phases and the quantum phase 
transitions among them. In this section, we present a clear evidence from the entanglement 
spectrum that the ground state of the SU(4) Heisenberg model (78) is indeed in an SPT phase 
protected by SU(4) [projective unitary group PSU(4), precisely^*] symmetry. 


3.3.1. Haldane phase -the simplest example 

To understand the structure of topological phases in the case of SU(7V) symmetry, it is con¬ 
venient to begin with the simplest prototypical case N = 2. Since the Haldane’s conjecture, we 
know that the ground-state properties of the spin-S Heisenberg chain are qualitatively different 
depending on the parity of 2S [118, 1 19]; when 2S = even, the ground state is in a featureless 
non-magnetic phase {Haldane phase) with the gapped triplon (S — 1) excitations in the bulk, 
while, for odd 25', we have a gapless (i.e., algebraic) ground state with spinon (S = 1/2) exci¬ 
tations. This conjecture has been later confirmed both by the construction of a rigorous example 
(see below) [140, 141] and by extensive numerical simulations [142-144]. Soon after, it has 
been pointed out that the featureless gapped ground state of the integer-S' spin chains may have a 
hidden “topological” order characterized by non-local order parameters [145-148] at least when 
S is an odd integer [149]. However, it was not until the concept of SPT phases was established 
that the true meaning of “topological order” in the Haldane phase was fully understood [37]. 
Now it is realized that the gapped phases in integer-spin chains with some protecting symmetry 
(e.g., time-reversal, reflection) are further categorized into topological phases and trivial ones. 

To understand the properties of the Haldane phase of integer-spin antiferromagnets, it is 
convenient to consider the spin-1 VBS state introduced by Affleck, Kennedy, Lieb, and Tasaki 
[140, 141]. The basic idea of construction is to first decompose an S' = 1 at each site into 
a pair of S = l/2s, then form uniform tiling of dimer singlets (hence ‘valence-bond solid’) 
among S=l/2s on neighboring sites, and fuse the S = 1/2 pairs at the same site back to the 
original spin-ls (see Fig. 9). One of the most convenient ways of representing the VBS state is 
the matrix-product-state (MPS) representation [150-152]'^: 


|S = lVBS)„,^ = (g)A 

i 

= X! {A{mi) ■ ■ ■ A{'mi) ■ ■ ■ A{'mL)}a,p\'mi) ®\'m2) ® ■ ■ ■\rnL) , 

mi—x,y,z 


( 88 ) 


where {| a;), | y), | z)} are defined by 

1 ^). = -^(1 + 1 ).- I - 1)0 

|y), = ^(| + l), + |-l),) 
\z)i = |0)i 


(89a) 

(89b) 

(89c) 


^^The difference between SIJ{N) and PSU(A^) may become clear below. 

^^The symbol (3) j implies both the matrix multiplication and the tensor product of the local (spin-1) Hilbert spaces. 
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and the matrices Ai and A(a) (a = x, y, z) are given in terms of the Pauli matrices; 


A{a) = (t“ (a = X, y, z) 

Ai= ^ A{a)\a)i = a'^\x)^ + ay\y)i+a^\z)^ . 

a—x,y,z 


( 90 ) 


Remarkably, it can be shown [140, 141] that this quantum many-body state (88) is the unique 
ground state^° of the following simple spin-1 Hamiltonian (VBS model): 

= E (s* • S*+1 + ^ (S* ■ • (91) 

It is easy to see that the VBS state (88) describes a non-magnetic short-range phase with 
an excitation gap. To see whether the state is topologically non-trivial or not, it is useful to 
consider how the state (88) transforms under the symmetry operation. Being non-magnetic, the 
bulk does not respond to the symmetry operation but the edges do. As the consequence, the 
symmetry operation gets fractionalized into two pieces; one acts on the left edge and the other 
on the right. For instance, the ground state of the spin-1 AKLT model (91) 15 = 1 VBS)q, 
hosts two emergent 5 = | spins (i.e., a, jd =t,4.) on both edges and hence transforms under the 
SO(3) rotation as 

|5 = 1 VBS)„,;3 ^ = 1 VBS)„,,;3. , (92) 

ol' ,/3' 

where U is the S = ^ rotation matrix of SU(2). Putting it another way, U serves as the mathe¬ 
matical labeling of the physical edge states. It is important to note that U appearing in Eq. (92), 
in general, may be a projective representation of SO(3) as both W and U appear in the equation. 
Since in the VBS state (88) U belongs to a non-trivial projective representation that is intrinsi¬ 
cally different from any irreducible (integer-spin) representations of the original SO(3), one sees 
that |5 = 1 VBS)q,^^ is in a non-trivial topological phase with emergent edge states. 

On the other hand, one can construct another model state of a spin-1 chain; 

|5 = lVBS-II)„,;3 = (g)Sz (93) 

i 

with the following 3x3 matrices Be 

B, = S^\x), + Sy\y), + S^\z), 

/O 0 0 \ / 0 0 A 

5 ^ = 0 0 , 5 ^ = 0 0 0 , 5 ^ 

\0 i 0 / 0 0/ 

By construction, it is obvious that the above state also exhibits edge states with spin-1, and one 
may suspect that it describes a SPT state. Using [5“]bc = —one can easily see that the state 
(93) transforms as before [see Eq. (92)] but with U now belonging to the spin-1 representation. 



(94) 


^**11 is unique in a periodic system or in the thermodynamic limit. On a finite open system, the ground state is four-fold 
degenerate because of the edge modes [141]. 
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Since the spin-1 representation is trivial in the sense of projective representation of SO(3), one 
can eliminate the would-be edge states by continuously deforming the Hamiltonian [134] and this 
ground state is in a trivial phase. This reasoning may be readily generalized; when U transforms 
like a half-odd-integer spin, the phase is topological, while when U transforms in an integer-spin 
representation [i.e., linear representation of SO(3)], the system is in a trivial phase. 

For later convenience, we rephrase the situation in terms of Young diagrams. The spin-^ 
representation of SU(2) is represented by the following Young diagram made of 2S boxes; 

m-'-D- (95) 

V V 

S/' 

2S boxes 

Then, the above result may be summarized as follows; when U belongs to the representations 

□ (96) 

the state represented by the corresponding MPS is topologically non-trivial as we cannot annihi¬ 
late these “emergent” edge spins by fusing physical integer spins on neighboring sites. 

On the other hand, the phase is trivial when 


m , I I I 11 ,--- • (97) 

That is, the number of boxes (mod 2) in the Young diagram for the representation to which U 
belongs labels the topological classes protected by SO(3) and leads to the Z 2 classification of the 
SO(3) SPT phases [153, 154]^'. As will be seen in Sec. 3.4, the N = 2 models (both g-e and 
p-band) host several different Haldane phases (in the spin, orbital, and charge sectors) in their 
phase diagrams. 


singlet pair 




m □□ 

spin-1 spin-1 


Figure 9: (Color online) Valence-bond construction of the spin-1 VBS states (88). Two S = 1/2 spins (ancillary 
qubits) at each site are symmetrized to obtain physical spin-1. If we replace the spin-l/2s with spin-ls and apply the 
anti-symmetrization at each site, we obtain the state (93). 


In the argument presented here, the Haldane phase is protected by the on-site symmetry SO(3) [not SU(2)]. However, 
it is known that other discrete symmetries, e.g., time-reversal and Z 2 XZ 2 can also protect the Haldane phase [133, 134]. 
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3.3.2. SU(N) topological phases 

Using the MPS representation [155] of the gapped ground state in one dimension, the above 
“physical” idea can be generalized and made mathematically precise. In fact, when a given 
ground state that is represented by an MPS 

A{mi)A(m2) ■ ■ ■ A(mL)|mi)(g) • • • <^\mL) (98) 

{mi} 

is invariant under some symmetry G, a 79-dimensional unitary matrix Ug (g G G) exists such 
that [156] 

A(m,) , (99) 

where A{mi) denotes the 79x79 MPS matrices corresponding to the local state \mi) and is 
a phase that depends on g G G. As has been mentioned above, the unitary matrix Ug is in fact a 
projective representation of the symmetry G, that corresponds to the physical edge states [133]. 
Therefore, the enumeration of topologically stable phases in the presence of symmetry G boils 
down to counting the possible (non-trivial) projective representations of G [153]. This problem 
was solved for SU(7V) and other Lie groups in Ref. [122] and the picture in the previous section 
basically generalizes to the case of SU(7V) with some mathematical complications. 

Now the role of SO(3) in the previous section is played by PSU(A^) ~ SU(7V)/Zjv [note 
SO(3) ~ PSU(2)]. Considering PSU(7V), instead of SU(7V), amounts to restricting ourselves 
only to the irreducible representations of SU(7V) specified by Young diagrams with the number 
of boxes ny divisible by N [i.e., ny = Nk (A: = 0,1,.. .)]^^. This subset of irreducible represen¬ 
tations roughly corresponds to the integer-spin ones in the SU(2) case. In view of the results of 
the Lieb-Schultz-Mattis argument presented in Sec 3.2.4, considering the symmetry PSU(7V) is 
quite reasonable for the hunt for gapped topological phases in one dimension, where no genuine 
topological phase with topologically degenerate ground states exists [39]. 

As in the previous section, the topological class of a given ground state (typically written as 
an MPS) is determined by looking at to which projective representation the unitary Ug of the state 
belongs. Since inequivalent projective representations of PSU(A^) are labeled by ny (mod TV) 
[122], there are TV — 1 non-trivial topological classes specified by the Z^r label ntop = ny (mod 
TV). In the following, we use the name “class-ritop” for the topological phase with ntop = ny (the 
class-0 corresponds to trivial phases). 

We can think of other protecting symmetries associated with SU(TV) [122]. For instance, 
when we take SU(4)/Z2 ~ SO(6), we consider only the SU(4) irreducible representations with 
ny = 0 (mod 2), which may be viewed as the linear representations of SO(6) (containing no 
spinor representation). In this case, we are led to a Z 2 -classification. 

A remark is in order about the definition of the topological class. In contrast to the SU(2) 
case where all the irreducible representations are self-conjugate, we must distinguish between an 
irreducible representation and its conjugate in SU(TV) (see Fig. 10). The relation (99) suggests 
that if we have the edge state transforming under the projective representation 7Z [of PSU(TV)] 
on the right edge, we necessarily have its conjugate 7Z on the left. This means that when we 
talk about the topological class we must first fix which edge state we use to label the topological 
phases. Throughout this paper, we define the topological class by the right edge state [i.e., by Ug 
acting from the right in Eq. (99)]. We will see, in the next section, that the SU(TV) VBS states to 
be discussed in Sec. 3.3.3 belongs to class-TV/2. 


refer the readers who want to know more about the mathematical details to Sections 11 and III of Ref. [122], 
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conjugate 



(a) 



Figure 10: (a) Young diagrams for an SU(Af) representation and its conjugate for N = 4. (b) Examples of self-conjugate 
representations {N = 4). 


3.3.3. VBS construction of the model wave functions 

In the previous section, we have seen that there are N topologically distinct phases in the 
presence of PSU(-/V)(~ SU(A^)/Zjv) symmetry [122]. To understand their physical properties, 
it is convenient to analyze the corresponding ‘hxed-point’ states that are easy to analyze. In the 
Haldane phase of integer-spin antiferromagnets, the spin-1 VBS state (88) discussed in Sec. 3.3.1 
does the job. 

A similar strategy applies to SU(A^) {N\ even, N > A.\ see Fig. 11) to construct the SU(A^) 
generalization of the VBS state. For instance, we use a pair of six-dimensional representations 
(Q) to build the SU(4) VBS state shown in Fig. 1 1(a) [103, 157]. At the last stage of the construc¬ 
tion, we project the tensor product 000 onto the 20-dimensional representation to obtain the 
physical Hilbert space. The resulting state may be conveniently expressed in the form of the 
MPS (98) using 6x6 matrices A(m) [158]. The parent Hamiltonian of this SU(4) VBS state 
reads as^^ [103, 157] 

^vfif = E {‘5* • '^*+1 + ^ > (100) 

where we have introduced a short-hand notation Si ■ Similarly, from 

a pair of 20-dimensional representations ^ of SU(6), we can construct the VBS ground state 
of the following Hamiltonian [with Si belonging to the 175-dimensional representation |^; see 
Fig. 11(b)] [158]24 


n 


N=6 

VBS 


E|'^* • ‘Si+i -f 

i ^ 



{Si ■ + 


1 

18288 


{Si ■ 



( 101 ) 


To be specihc, let us restrict ourselves to the case of SU(4) to describe the main features of 
the VBS wave function. Most of the properties are carried over to general SU(A^) (N: even) with 


fact, the parent Hamiltonian contains two free parameters aside from the overall factor. Requiring that 
(Si • Si^i)‘^ and (Si • should not appear, we obtain Eq. (100). 

After completing this paper, we were informed that T. Quella et al. had independently obtained similar results in 
their unpublished work [159]. 
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Figure 11: (Color online) Valence-bond constnaction of the topological VBS states, (a) and (a)’: class-2 SU(4) VBS 
state, (b) and (b)’: class-3 SU(6) VBS state. From Ref. [103]. 


due modifications. The ground state is SU(4)-symmetric and featureless in the bulk, and has the 
‘spin-spin’ correlation functions [103] 


( 5 , 




j+n 



n ^ 0 
n = 0 


( 102 ) 


that are exponentially decaying with a very short correlation length 1/ In 5 « 0.6213. Neverthe¬ 
less, the system hosts emergent edge states near the boundaries. In fact, if one measures (5/') 
(with 5/^ being any three commuting generators of SU(4), or equivalently, independent linear 
combinations of local fermion densities ^ iCga,i + ^Cea,i), one can clearly see the 

structure localized around the two edges (Fig. 12. See also the bottom panel of Fig. 19 for a 
similar plot for the original fermionic model). At each edge, there are six different states (i.e., 
0) distinguished by the value of the set of the three generators {Sf). This implies that Ug in 
Eq. (99) transforms like |^, telling us that the ground state of the VBS Hamiltonian (100) falls 
into the class-2 topological phase protected by the on-site PSU(4) symmetry. 

Corresponding to the physical edge states, the entanglement spectrum exhibits a particular 
degeneracy structure. In fact, the entanglement spectrum of the SU(4) VBS state consists of 
a single six-fold degenerate level, which leads to the von Neumann (bipartite) entanglement 
entropy S'vN = In 6.^^ This extremely simple structure is peculiar to the solvable VBS model 
(100) and, for a generic SU(4) Hamiltonian in the same class-2 topological phase, we expect 
the degeneracy structure compatible with the SU(4) irreducible representations having ny = 2 
boxes. Fig. 13 is the plot of the entanglement spectrum of the pure SU(4) Heisenberg model (79) 
obtained [158] by infinite time-evolving block decimation (iTEBD) method [160, 161]. One can 
clearly see that the level structure is totally consistent with the topological class-2. By linearly 
interpolating between the VBS model (100) and the pure Heisenberg model (79), we can show 
[158] that the overall structure of the entanglement spectrum is maintained indicating that both 
models belong to the same unique SPT phase. 


construction, it is obvious that for general N (even), the single entanglement level of the VBS state (Fig. 1 1) is 
W!/[(Af/2)!]^-fold degenerate. 
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Figure 12: (Color online) Local fermion density ria^i = i^ga,i + i^ea,i calculated for one of the 36-fold 
degenerate SU(4) VBS states on a finite open chain with 30 sites. There are six edge states as two of the four take 
the same value (and so do the other two) around the edge. Note that the local constraint ria.i = 4 is satisfied at 

each site. 


Before concluding the discussion of the SU(iV) SPT phases, a few remarks are in order. One 
is about the other topological phases. The group-cohomology classification tells us that there 
are two more topological classes in PSU(4)-invariant systems [122]. In fact, the VBS states 
studied in Refs. [141, 162-164] correspond to these two. The entanglement spectra calculated 
for these VBS states agree with those expected from the class-1 and 3 topological phases. Quite 
recently, the phases of SU(3)-invariant spin chains, including the SU(3)-counterpart of the above 
two phases were investigated in Ref. [165]. 

Also interesting is the characterization of these topological phases with non-local order pa¬ 
rameters. It has been known that the Haldane phase of SU(2) spin chains is characterized by a 
pair of non-local order parameters (string order parameters) [145, 147-149, 166]. These non¬ 
local order parameters are not only useful as ‘working indicators’ of the topological Haldane 
phase but also have intimate connection to the modern characterization of the ID SPT phases in 
terms of the projective representation Ug [see Eq. (99)] [167, 168]. 

This concept can be generalized to other symmetries [169-171]. For instance, a set of 
2(^N — 1) non-local ‘string’ order parameters distinguish between the N topologically distinct 
phases predicted by group cohomology [158]. The characterization of the SPT phases by non¬ 
local correlation functions seem quite interesting in view of the recent development in real-space 
imaging techniques [172] as some of the string order parameters are written only in terms of 
fermion densities which are, in principle, detectable in experiments. 

Last, we comment on the relation between the trap potentials and the edge states character¬ 
istic of SPT phases. In the usual setting of harmonic potentials, it is known [173, 174] that an 
island of Mott insulating region (Mott core), that is surrounded by a compressible metallic region, 
is formed around the center of the trap; the above argument for SPT phases holds only within the 
Mott core. Due to the interaction between the Mott core and the metallic region surrounding it, 
the structures around the (smeared) edges may be substantially reduced^® [106, 107]. However, 


^^Especially in the case of, e.g., the charge Haldane phase (see Sec. 3.4.1 and Fig. 16 for the definition) where the 


39 






14 


12 

10 



50 



Q t.^. 

0 20 40 60 80 100 120 

a 


Figure 13: (Color online) Entanglement spectrum of the SU(4) Heisenberg model (78). From Ref. [158] 


thanks to the recent success in creating a box trap [175], one may resolve this problem. Com¬ 
bining the box trap and the technology of single-site detection (see [176] for a recent review), 
quantum simulation of various SPT phases [including our SU(A^) ones] will become feasible in 
the near future. 

3.4. Example of phase diagrams 

We will now turn to the presentation of some numerical phase diagrams that were obtained 
using DMRG simulations on finite chains of length L with open boundary conditions. Since 
we focus on particular phases occurring at fixed densities, we do not consider the (harmonic) 
trapping potential. In real experiments, detection could be achieved using local quantities. Let 
us emphasize that we will discuss mostly gapped phases, so that they are stable with respect 
to small perturbations and should be more robust for observation. As to the simulations, let us 
simply mention that we have only implemented the U(l) quantum numbers corresponding to 
the separate conservation of the fermion numbers of each color, though it should be possible to 
ease simulations by implementing the full SU(7V) symmetry of the models [177]. We refer to 
Ref. [103] for more technical details. 

For the sake of clarity, in this section, we restrict ourselves to the half-filled cases which 
already exhibit a rich variety of phases, including topological ones. Away from half-filling, one 
expects on general grounds the occurrence of gapless Luttinger liquid-like phases with (diagonal) 
density or superconducting (i.e. off-diagonal) correlations [178] or degenerate Mott phases [179] 
for some other (in)commensurate fillings. 

For simplicity, when considering the g-e model (34), we will assume that both orbitals are 
equivalent, so that they have identical hoppings tg = = t, chemical potentials pg = pe = 

and interactions Ugg = Uee = Umm- Therefore, we will consider neither the case with spin- 
imbalance where spin-polarization effects may dominate [107] and give rise to FFLO physics. 


edge states directly couple to the trapping potential, the inhomogeneity of the potential causes destructive effects. 
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nor the strongly anisotropic case where one of the orbitals would be much more localized than 
the other (e.g., tg <C tg) hence giving rise to physics of the SU(A^) Kondo lattice model [17, 180- 
184], 


3.4.1. N = 2 g-e model 

We present in Fig. 14 typical phase diagrams of the g-e model (34) with N = 2 which exhibit 
a large variety of phases: (i) charge density wave (CDW), (ii) orbital density wave (ODW), 
(iii) spin-Peierls (SP), (iv) charge Haldane (CH), (v) orbital Haldane (OH), (vi) spin Haldane 
(SH), and (vii) rung singlet (RS) (see Fig. 16 and Table 3). The CH phase is the collective 
Haldane state formed by the (spin-singlet) charge pseudo-spin-1 states: (^gi(^g 2 ^ei^e 2 \^) = 

4), CgiCg2|0) + cL^LlO) (n = 2), and |0) (n = 0) [see Eq. (69b) for the definition of the 
charge pseudo-spin] and is characterized by the existence of non-local string order in the charge 
distribution. Clearly, it requires strong charge fluctuations and is not a Mott insulator. Similar 
‘Haldane states’ in the charge sector have been found in the study of related multi-component 
fermions [45, 51, 107]. As CH is a collective insulating state with a charge gap analogous to 
the spin gap in the Haldane-gap systems, it, together with its bosonic counterpart, is also dubbed 
Haldane insulator in literatures [185-187]. To give a simple picture of these phases, we provide, 
respectively in Figs. 15 and 16, a cartoon of the two-fold degenerate possible density waves 
(CDW and ODW) and the non-degenerate Mott phases (SH,CH, and RS) that we have found. 
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Figure 14: (Color online) Phase diagram of the g-e model with N = 2 at half-filling. Left and right panels correspond 
respectively to = —1 and t4f '’/f = 1, which allows to realize various phases: spin-Peierls (SP), charge/orbital 

density waves (CDW/ODW), rung singlet (RS), charge-Haldane (CH), orbital-Haldane (OH) and spin-Haldane (SH). 


These very rich phase diagrams are in rather good agreement with the low-energy predictions, 
and they were already discussed in Refs. [103, 105]. In Fig. 14, one notices that the phases 
concerning the charge sector (CDW and CH) and those concerning the orbital sector (ODW and 
OH) appear in a very symmetric manner. In fact, this is a direct consequence of the following 
symmetry (a generalization of the Shiba duality) that the N = 2 g-e model possesses [103]: 


y ^ -y -f 


yg-< 


TZ9-e TT 
K pY 1 rr 


u„ 


( 103 ) 


that swaps a phase related to charge and the corresponding orbital phase. 
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Figure 15: (Color online) Two density-wave states for N = 2. In-phase and out-of-phase combinations of two density 
waves in g and e orbitals respectively form (a) CDW and (b) ODW. Adapted from Ref. [103]. 
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Figure 16: (Color online) Four translationally invariant insulating states for N = 2. (a) spin Haldane (SH), (b) orbital 
Haldane (OH), (c) charge Haldane (CH), and (d) rung-singlet (RS) phases. Singlet bonds formed between spins (orbital 
pseudo-spins) are shown by thick solid (dashed) lines (singlet bonds are not shown in (c)). Charge pseudo-spin is defined 
as Kf = ni/2 — 1 and the = 0 state is a linear superposition of and (state with = 0). 

Dashed ovals (rectangles) denote spin-singlets (triplets). Adapted from Ref. [103]. 


3.4.2. N = 2 p-band model 

We present here some data for the p-band model (63) with iV = 2 at half-filling. We refer the 
readers to Ref. [103] for the full phase diagram. We focus here along the line Ui = 3U2 which 
could be realized using the harmonic trapping [precisely, traps that are axially symmetric with 
respect to the chain direction; see Eq.(60)] and we simply plot some local quantities, that would 
be most easily measured experimentally using the existing quantum-optical techniques, namely: 
local densities and kinetic energy 


na{i) 





(104) 


In Fig. 17, one can clearly identify the edge states both for repulsive (Ui = SU 2 > 0) and 
attractive (Ui = 3 U 2 < 0) interactions. They correspond respectively to the induced edge states 
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local quantities 


in the spin and the charge Haldane phases (SH and CH; see Fig. 16) and are simple manifestations 
of the topological properties in the bulk (see the discussion in Sec. 3.3). The existence of the CH 
phase in fermionic systems was pointed out in Refs. [45, 51] for a related but different multi- 
component Fermi system. The CH phase in the p-band model was first found in Ref [107]. 

Let us remind the readers that, for this model, the sign of the interaction can be flipped 
formally [(C/i, C/ 2 ) —> {—Ui, —C/ 2 )] using spin-charge interchange transformation [107] similar 
to the one used in Sec. 2.3.3 to discuss the attractive (U < 0) side of the half-filled single¬ 
band SU(2) Hubbard chain (see also the discussion in Sec. 3.1.2). Because of this spin-charge 
symmetry in the N — 2 p-band model, SH and CH appear in a symmetric manner in the phase 
diagram (see Fig. 13 of Ref. [103]). 



J^_L 

64 96 

i 




Figure 17: (Color online) Local quantities (fermion density and kinetic energy) for the p-band model at half-filling with 
N = 2 obtained on L = 128 chain using DMRG. Left panel: (17i/f, (72/f) = (12,4) exhibits spin edge states 
characteristics of SH phase. Right panel: (Fi/f, C/ 2 /t) = (—12, —4) exhibits charge edge states characteristics of CH 
phase. 


3.4.3. N = 4 g-e model 

We now turn to the = 4 g-e model (34) and present its phase diagrams for several fixed 
values of jt in Fig. 18. Due to the existence of the symmetry V —>■ —V that swaps charge 
and orbital when = 0 [103], the phase diagram shown in the second panel is symmetric 
with respect to C = 0. 

In the weak-coupling regime, we only find conventional degenerate phases with broken trans¬ 
lation symmetry (CDW, ODW and SP; See Table 3 for the meanings of the abbreviations), that 
are in good agreement with the low-energy predictions [103]. However, in the intermediate and 
strong-coupling regions, we find the SU(4) topological phase as has been predicted by the strong¬ 
coupling argument (see Secs. 3.2.2 and 3.3.3). Its signatures are again given by the existence of 
non-trivial edge states [see Fig. 15(a) of Ref [103], which looks similar to Fig. 19(c)] or could 
also be probed by computing the degeneracy of the entanglement spectrum [158] [see Fig. 13 
for the entanglement spectrum of the effective model (78) that describes the topological phase]. 
There is thus a quantum phase transition between (weak-coupling) SP and the SU(4) topological 
phase (class-2 SPT phase) in the strong-coupling region. The nature of the universality class of 
the transition is, however, difficult to determine numerically. In Ref. [157], it was conjectured 
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that the quantum phase transition is governed by an SU(4)2 CFT with the central charge c = 5. 
The high value of the central charge calls for large-scale numerical simulations. 
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Figure 18: (Color online) Phase diagram of the g-e model with = 4 at half-filling. From top-left to bottom, we 
have fixed 0 and 1, which allows to realize conventional phases (SP, CDW and ODW) as well as the 

topological SU(4) one. The symmetry with respect to V = 0 in the second plot is a direct consequence of the orbital- 
charge interchange that maps CDW into ODW. See Table 3 for the meaning of the abbreviations ‘ODW’, ‘RS’, etc. 


3.4.4. N = 4 p-band model 

When considering the TV = 4 p-band model, we restrict ourselves to the case with half¬ 
filling. Here we only plot relevant local quantities (local fermion densities and kinetic-energy 
density) in Fig. 19 for several couplings along the line Ui = 3(72 (harmonic potential). The full 
U 1 -U 2 phase diagram obtained in Ref. [103] is also shown in Fig. 20. The definitions of the local 
quantities are similar to Eq. (104) with the only difference that a = 1,..., 4 for = 4. At weak- 
coupling, Figs. 19(a,b) indicate the existence of CDW phase for attractive interactions {Ui < 0) 
and SP in the repulsive (Ui > 0) case, as predicted using the low-energy field theories [103]. For 
strong repulsive interactions. Fig. 19(c) shows a very different behavior since non-trivial edge 
states appear while the bulk becomes featureless, and this is in agreement with the property of 
the proposed topological SU(4) phase (see Fig. 12). Indeed, each edge state is 6-fold degenerate 
as it results from choosing 2 colors among 4, and this degeneracy can be directly seen in the 
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entanglement spectram (see the lowest level in Fig. 13). As in the iV = 4 g-e model, we expect 
a quantum phase transition in the SU(4)2 universality class between SP and SU(4) topological 
phases. 


phases 

abbreviation 

SU(V) 

orbital (T) 

spin-Haldane 

SH 

S'= 1 

local singlet 

orbital-Haldane 

OH 

local singlet 

N/2 

charge-Haldane 

CH 

local singlet 

- 

orbital large-D^; j, 

OLDa: y 

local singlet 

N/2 

rung-singlet (OLD^) 

RS 

local singlet 

N/2 

spin-Peierls 

SP 

- 

N/2 

charge-density wave 

CDW 

local singlet 

local singlet 

orbital-density wave 

ODW 

local singlet 

N/2 


Table 3: List of dominant phases and their abbreviations. Local SU(A^)/orbital degrees of freedom are shown, too. 
(see also Fig. 16) ‘SH’ and ‘OH’ appear only in the N = 2 case. ‘RS’ and ‘ODW’ respectively are a ‘large-D’ state 
(g)|T^ = 0) and the ‘Neel-ordered’ state of T = N/2 chains. 
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local quantities 



Figure 19: (Color online) Local quantities for the p-band model at half-filling with = 4 obtained on L = 64 
chain using DMRG. The three panels correspond to several points relevant for harmonic trapping: (a) (C/i/f, U 2 lt) = 
(-3,-1) corresponding to a CDW phase; (b) {Ui/t, C/ 2 /f) = (3,1) exhibits strong dimerization, as expected for a SP 
phase; (c) {Ui/t,U 2 /t) = (9,3) shows non-trivial edge states characteristics of the topological SU(4) phase. Note that 
in the first two panels, we only plot data in the bulk of the chain for readability. 
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Figure 20: (Color online) Phase diagram for half-filled N = A p-band model (63) obtained by DMRG on L = 32. 
Dashed line corresponds to the condition Ui = 3(72 satisfied for an axially symmetric trap. U 2 = 0 correspond to two 
decoupled SU(4) Hubbard. Adapted from Ref. [103]. 
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4. Concluding remarks 

Alkaline-earth and ytterbium cold atomic gases allow the realization of fermionic cold gases 
with SU(A^) symmetric interactions in a very controlled way. On top of the almost perfect SU(iV) 
symmetry, one has the advantage of being able to use two orbital levels so that we can play 
with spin-orbital coherent exchange interactions, as demonstrated experimentally [26-28]. Both 
ingredients are expected to lead to a variety of states of matter and phenomena depending on the 
microscopic parameters, the dimensionality and so on [17]. In such a context, our review has 
focused on the one-dimensional case, where powerful analytical and numerical techniques are 
available. 

Our review has covered, in depth, the case of single-band SU(A^) Fermi-Hubbard model 
which ultracold alkaline-earth fermions loaded into a one-dimensional optical lattice can simu¬ 
late up to N as large as 10. The physics of the single-band Fermi-Hubbard model with a single 
contact interaction U is very rich for general hllings and turns out to be very different between 
N = 2 and N > 3. One fascinating aspect of the model is the possible occurrence of the (hnite- 
U) Mott transition for commensurate hllings when > 2 in sharp contrast to the N — 2 case 
where the critical value is [4 = 0 and additional nearest-neighbor interactions are necessary to 
have a hnite [4- For one atom per site, a hlling at which we have least three-body losses, the low- 
energy physics of the Mott-insulating phase is described by the integrable Sutherland model with 
— 1 gapless relativistic modes. For incommensurate hllings and in the low-density regime, 
a new Luther-Emery liquid with a gap in the SU(Al) sector and one gapless bosonic (charge) 
mode emerges in the attractive U < 0 case. The instability toward a paired fermionic super- 
Huid (an analogue of the BCS superconductivity) is completely suppressed in this phase and the 
dominant superhuid instability occurs in the A^-particle channel, e.g., the trionic and quartetting 
instabilities when N = 3 and A^ = 4, respectively. Finally, N > 2 systems support fully gapped 
Mott-insulating phases for special commensurate hllings; these phases display a bond-ordering 
with ground-state degeneracies as in a SP phase at half hlling. 

All these results obtained for the single-band case clearly imply a rich variety of physics when 
N > 2 which cannot be realized in the SU(2) two-component Fermi gas with a contact inter¬ 
action. However, as already seen, all the fully gapped SU(A^) Mott-insulating phases occurring 
in the single-band model spontaneously break the translation symmetry. In this respect, many 
interesting fermionic or bosonic SPT phases are beyond the scope of the single-band SU(A^) 
Fermi-Hubbard model. The realization of these important phases thus necessitates the general¬ 
ization of the model by introducing, for instance, additional degrees of freedom. To this end, 
we have presented two microscopic models which describe the low-energy physics of the short- 
range interacting two-orbital SU(A^) cold fermions on a lattice, namely the g-e model and the 
p-band one. Both models can be realized with ultracold alkaline-earth or ytterbium atoms. The 
zero-temperature phase diagrams of these models are very rich and harbor very different Mott 
and collective insulating phases. On the basis of low-energy held-theory approaches as well as 
rigorous constructions of the ground states in some strong-coupling regime, a variety of phases 
have been found that are either degenerate (ODW or CDW, SP etc.) or non-degenerate (triv¬ 
ial singlet, various Haldane-like phases, as well as other SPT phases). All these results can be 
conhrmed numerically by mapping out various phase diagrams, as we have presented for even 
N = 2 and 4 at half-hlling for both models. 

Let us emphasize that in the SPT phases, the edge states are protected (and thus cannot be 
removed without closing a gap) as long as a given symmetry is present. This is precisely the 
case, for instance, for the SU(Af)-protected topological phase that we have discussed here and 
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our observation provides one microscopic realization of one of the — 1 possible SPT phases 
predicted in Ref. [122]. Therefore, its experimental observation using, for instance, polarization 
measurements [188, 189], site-resolved imaging technologies [172], and spin-selective detec¬ 
tion [190] looks like an exciting possibility, although the edge states may be suppressed or even 
absent if one takes into account a harmonic trap [106, 107]. In order to circumvent this diffi¬ 
culty, an interesting framework would be to use a box-shaped trapping potential [175] where 
presumably edge states should be more visible. 

Of course, there are many important topics which have not been covered in this review. For 
instance, when discussing the two-orbital models, we have assumed that both orbitals are not very 
different from each other, i.e., atoms have similar hoppings and chemical potentials for both g 
and e orbitals, as can be achieved experimentally using, for instance, the p-band levels discussed 
in Sec. 3.1.2. Therefore, we have not touched upon, e.g., the fascinating Kondo physics in the 
context of ultracold alkaline-earth atoms [181-184]. Also, we have not discussed interesting 
phases that would appear away from half-hlling. For instance, as in the single-band case, we 
may expect various superfluid instabilities to occur for incommensurate fillings. This problem is 
yet to be fully understood except for the N = 2 case where we found several competing pairing 
phases as well as Mott ones [178]. 

Recently, it has also been proposed that using synthetic gauge field on one-dimensional 
SU(A^) cold atoms with contact interactions may be related to 2D Chern insulators [191] or 
quantum Hall phases [192], which also paves the way to realize exotic phases of matter with 
these systems. In the light of the recent experimental achievements with alkaline-earth cold 
fermionic gases, we hope that it will be possible in the future to unveil part of the richness that 
we highlighted in this review. 
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